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A  COASTAL  DISPERSION  MODEL 
FOR  EFFLUENT  PLUMES 


SUMMARY 


A  computer  simulation  model   was  developed  to  plot  two  dimensional 
dilution  contours  of  an  outfall     discharging  into  the  time  variant 
coastal   currents  of  a  large  lake.     The  computer  model    is  based  on  an 
analytical    steady  state,  depth  integrated,  two  dimensional,  Gaussian 
cross-plume  model,  which  generates  successive  realizations  of  the 
farfield  plume.     The  current  speed  and  direction  data   (from  the 
current  meters)  are  used  to  map  dil ut ion /concentration  contours  in 
polar  coordinates.     The  model   assumes  that: 

-  the  effluent  field  moves  with  the  prevailing  currents  without 
disturbing  the  existing  flow  pattern  of  the  coastal   zone. 

-  the  diffusion  in  the  flow  direction  is  negligible  compared  to 
the  advection. 

-  the  lateral   eddy  diffusivity  is  a  function  of  the  initial    plume 
width. 

-  the  vertical   diffusion  is  negligible  and  the  effluent  is 
uniformly      distributed  over  the  available  water  depths. 

The  following  steps  produce  dilution  contours  from  a  diffuser 
outfall. 

-  determine  the  initial   dilution  due  to  the  jet  rise  using  the 
techniques  discussed  in  this  report; 

-  classify  the  current  episodes  according  to  the  speed  and  the 
direction  ranges; 

-  divide  each  episode  into  hourly  speed  and  direction  events; 
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-  calculate  the  farfield  dilution   (or  the  total   dilution   including 
the  initial   dilution)   for  each  hourly  event  at  all   points  of  a 
polar  grid; 

-  repeat  the  procedure  for  other  hourly  events  and  then  average 
the  results  over  the  whole  episode. 

These  successive  plumes  provide  statistical   estimates  of  the  plume 
transport  under  selected  current  regimes. 

The  model   was  applied  to  the  lakeview  Water  Pollution  Control    Plant 
(WPCP)  diffuser  outfall    in  Lake  Ontario,   using  the  1970-72  current 
meter  data,  to  predict  dilution  contours  of  a  conservative 
pollutant.     This  method  can  also  be  applied  to  non-conservative 
pollutants,  if  their  decay  rates  are  known.     The  model   can  be 
applied  to  any  diffuser  outfall    in  a  large  lake,  provided  it  is 
calibrated  by  collecting  the  current  meter  and  plume  tracking 
information,    simul tanously. 

Present  experience  suggests  two  improvements  to  the  computer  model. 

1.  As  the  model   was  found  to  be  more  sensitive  to  the  current 
direction  than  the  current  speed,  the  accuracy  of  diffusion 
field  calculations  can  be  improved  by  decreasing  the  current 
direction  segments  of  the  current  episodes. 

2.  Short  duration  current  episodes,  where  the  plume  length  is 
limited  by  advection  rather  than  diffusion,  should  be  included 
in  the  analysis. 
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A  COASTAL  DISPERSION  MODEL 
FOR  EFFLUENT  PLUMES 


1.   INTRODUCTION 

One  of  the  objectives  of  the  Great  Lakes  Water  Quality  Agreement 
(International  Joint  Commission  -  IJC)  of  1978,  signed  between 
Canada  and  the  United  States  of  America,  is  to  achieve  and  maintain 
mutually  acceptable  water  quality  conditions  in  the  Great  Lakes.   It 
was  recognized  at  the  time  that  it  might  not  be  economically  or 
physically  feasible  to  achieve  the  same  water  quality  in  the 
irmtediate  vicinity  of  an  effluent  outfall  as  in  the  remainder  of  a 
lake.  Consequently,  in  Annex  2  of  the  Agreement,  the  appropriate 
governing  bodies  were  given  the  responsibility,  with  limiting 
guidelines,  to  define  certain  areas  as  limited  use  zones  (LUZ). 
Within  the  LUZ,  the  water  quality  may  not  comply  with  the  water 
quality  objectives.  However,  before  any  water  passes  through  a  zone 
boundary,  it  must  achieve  the  designated  objective  of  the  receiving 
lake.  The  purpose  of  this  report,  which  is  a  joint  effort  by 
the  National  Water  Research  Institute  (NWRI)  and  the  Ontario 
Ministry  of  the  Environment  (M0E)}  is  to  describe  how  the  coastal 
currents  and  dispersion  characteristics  derived  from  the  existing 
climatological  data  bases  of  NWRI,  MOE,  and  Ontario  Hydro  can  assist 
in  delineating  the  LUZ  perimeter  of  an  outfall  discharging  into  a 
large  lake. 

Sewage  outfalls  and  other  discharge  outlets  are  usually  located  at 
lake  bottom  within  2  km  of  the  shoreline.  These  outfalls  are 
generally  multiport  diffusers  placed  such  that  the  diffuser  jets  are 
perpendicular  to  the  predominant  current  direction.  After  the 
initial  jet  mixing,  which  takes  place  within  a  few  minutes  or  a 
short  distance,  the  effluents  effectively  form  a  line  source  at  the 
surface.  Further  diffusion  of  this  line  source  of  effluents  is 
governed  by  the  prevailing  coastal  currents  and  their  eddies.  The 
dilution  will  depend  upon  the  location  of  the  effluent  plume  within 
the  coastal  zone,  current  transport  and  dispersion  characteristics 
of  the  area. 

The  coastal  currents  and  their  dispersion  characteristics  can  be 


obtained  by  direct  measurement  or  inferred  from  existing  data  bases 
after  site  specific  validation.     It   is  not  yet  possible  to  model   all 
of  the  complex  flows  that  have  been  observed  in  the  coastal    zones  of 
the  Great  Lakes  due  to  lack  of  knowledge  of  their  interactibns. 
Fortunately,  the  important  and  frequent  case  of  longshore  pollutant 
transport  and  dispersion  by  a  nearly  steady  shore  parallel   current, 
can  be  predicted  by  the  following  combined  statistical -analytic 
analysis.     Similar  shore  parallel    or  antiparallel  current  episodes 
are  selected   from  the  climatological   data  base.     (If  a  person  stands 
on  the  lakeshore  and  looks  towards  the  lake,   his  right  hand 
direction  is  defined  as   'shore  parallel'   and  his  left  hand  direction 
is  defined  as   'antiparallel').     These  are  used  as  the  input  data  to 
a  steady  state,  depth  integrated,  two-dimensional    Gaussian 
cross-plume  model  which  generates  successive  realizations  of  the 
effluent   (farfield)  plume.     The  minimum,  mean,  and  maximum  limits  of 
the  farfield  plume,  or  a  constant  dilution  contour  within  the  plume, 
provide  statistical1  estimates  of  the  plume  transport  and  dispersion 
under  the  selected  current   regime. 

A  computer  program  was  developed  to  carry  out  the  above  procedure 
for  conservative  and  nonconservative  substances.     This  program  of 
farfield  dilution  requires  a  value  of  nearfield  or  initial   dilution 
as  an  input,  which  can  be  calculated  separately  or  by  means  of  a 
user  subroutine.     Providing  the  climatological    data  base  and 
selected  current   regimes  are  representative  of  the  site  under 
investigation, a  reasonable  estimate  may  be  made  of  the  frequency, 
duration,   and  dimensions  of  future  effluent  plumes.     The  model    can 
also  be  applied  to  conditions  when  the  currents  are  in  directions 
other  than  shore  parallel    or  antiparallel,  for  example,  offshore  and 
onshore  current  directions. 


2.     THE  COASTAL  ZONE 

It  has  been  known  for  the  last  decade  that  at  the  edge  of  the  Great 
Lakes  there   is  a  5  -  10  km.  coastal   zone  where  the  lakewide 
circulation  adapts  to  the  presence  of  the  shore,   i.e.  the  lake 
currents  are  influenced  by  bottom  and  the  friction  shore  boundary 
(Csanady  1970,   1972a,   1972b).     Thus,  for  our  purposes,  the  coastal 
zone   is  defined  as  that  area  extending  offshore  over  which  bottom 
friction  effects  modify  the  lakewide  circulation  currents  to  the 
extent  that  the  predominant  local   currents  are  shore-parallel   and 
are  increasingly  dependent  on  local   wind  forcing   (Bull   and  Murthy, 
1978).     During  the   sunmer  months,  when  the  lakes  are  strongly 
stratified,  the  coastal    zone  is  narrower  than  during  the  winter 
months  of  weak  stratification   (Bull   and  Murthy,   1980).     Extensive 
analysis  of  NWRI  measurements  suggests  that  the  coastal    zone  can  be 
divided   into  two  distinct   boundary  layers,   an   inner  or  frictional 
boundary  layer   (FBL)  dominated  by  shore  and  bottom  friction  and  an 
outer  boundary  layer  (OBL)  of  adjustment  of  the  inertial   component 
of  midlake  circulation   (Murthy  and  Dunbar,   1977). 

The  above  division   is  most  easily  illustrated  at  a  given   site  by 
plotting  the  mean  shore  parallel   currents  versus  distance  from  shore 
(Fig.   1)   or  by  directly  examining  stick  vector  plots  of  the  observed 
currents   (Fig.   2).     The  plots  indicate  that  the  most  efficient 
dispersal   of  pollutants,   associated  with  the  highest  currents  and 
turbulence,  is  obtained  by  locating  the  diffuser  at  the  FBL-OBL 
interface  1-4  km.  offshore.     This   is  seldom  cost  effective  but   it 
should  be  kept  in  mind  that  as  the  diffuser  approaches  the  FBL-OBL 
interface  from  the  shore  both  the  nearfield  dilution,  which  usually 
depends  on  depth,  and  the  farfield  dilution,  which  depends  on  the 
current   speed   and   turbulence,   are  likely  to  improve,   all    other 
factors  being  equal . 
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3.     COASTAL  CLIMATOLOGY 

Field  measurement  programs  designed  to  obtain  time-series  records  of 
coastal   currents  and  temperatures  over  periods  of  many  months  give  a 
time  history  of  the  various  flow  regimes  present  in  the  coastal 
zone.     Such  records  show  an  extremely  complicated  flow  situation 
varying  by  season  and  from  place  to  place  that  is  often  difficult  to 
analyse.     As  yet,   nearshore  models  are  insufficiently  comprehensive 
to  be  able  to  predict  the  complex  array  of  flow  regimes  encountered 
in  this  area.     Hence  one   is   reduced  to  employing  a   statistical 
description  of  the  measured  flow  properties.     A  climatology  of  the 
coastal   zone  is  built  up  by  cataloguing  the  relative  frequency  of 
occurrence  and  average  duration  of  important  flow  regimes  in  an 
area.     These  are  then  used  to  estimate  the  flushing  capability  and 
dispersion  of  the  local   currents.     It  must  be  emphasized  that  some 
of  the  flow  regimes  which  can  be  identified  are  not  easily  analysed 
with  existing  nearshore  models,  either  because  of  the  inherent 
complexity  of  the  flow  regime  or  due  to  insufficient  information 
about  the  associated  boundary  conditions  and  forcing  functions. 
Engineers  and  planners  should  always  be  aware  that  present  dilution 
estimates  are  based  on  those  subsets  of  the  available  current 
measurements  which  are  sufficiently  well    behaved  to  be  analysed  with 
the  existing  nearshore  models.     Unusual    and  complex  current  episodes 
are  nearly  always  excluded. 

Over  the  past  decade  a  large  data  base  of  coastal   current 
measurements  has  been  built  up  by  researchers  from  NWRI,  MOE, 
Ontario  Hydro,  and  the  University  of  Waterloo   (Csanady  and  Pade, 
1970,   1971,   1972;   Kohli   and  Palmer,   1973;  Kohli,   1976a,  1976b;  Kohl i 
and   Farooqui  ,   1980;   Bull    and  Farooqui  ,  1976;  Jordan  and  Bull,   1977; 
Murthy  and  Dunbar,   1977;  Kaiser  and  Bull,   1978a,   1978b;  Bull   and 
Kaiser,   1980;   Bull,  Kerman ,  and   Lilly,    1980).     An  examination  of 
these  data  sets   indicates  that  strong  shore-parallel   currents 
dominate  the  organized  flow  patterns  within  the  coastal    zones  of  the 
Great  Lakes.     For  the  present  calculations  the  current  climatology 
is  divided  into  two  sets.     The  first  contains  all   the  quasi -steady 
shore  parallel    flow  regimes  persisting  for  a  day  or  longer  which  are 
responsible  for  longshore  pollutant  transport. 


The  second  or  miscellaneous  set  contains  the  short  term,   irregular, 
and   random  flows  which  merely  distribute  the  effluent   in  the 
inmediate  vicinity  of  an  outfall.     This  division  is  necessary 
because  the  Gaussian  plume  model    is  only  valid   in  quasi -stationary 
current  regimes. 

The  longshore  current  flows  are  further  subdivided  into  two 
directions  (shore  parallel   and  antiparallel )   and  three  speed  ranges 
(0-5,   5-15,   >15  cm/s).     The  statistical   behaviour  of  the  effluent 
plume  can  now  be  modelled  by  calculating  the  farfield  concentration 
occurring   under  these  six  flow  regimes.     The  development  of  a 
persistent   longshore  current   regime   in  any  given  direction 
invariably  suggests  a  current  reversal.     However,  there  is  usually  a 
mass  exchange  between  the  coastal   and  offshore  waters  during  current 
reversals  (Csanady  1970).     Hence  the  plume  develops  in  clean  water 
and  the  preceding  flow  field  history  is  not   required.     Strong  shore 
parallel    currents  usually  persist  from  1-7  days  and  follow  the 
coastline.     At  an  average  speed  of  10  cm/s  a  conservative  effluent 
can  easily  be  carried  25  km.  or  more  from  the  source. 

During  the  winter  months  the  coastal   currents  exhibit  a  classic 
bi -modal    character,  flowing  either  parallel    or  antiparallel    to  the 
shore  (Fig.   2).     There  may  also  be  some  asymnetry  between  the  shore 
parallel    and  antiparallel   episodes,  the  currents  being  stronger  and 
of  longer  duration   in  one  direction  than  the  other  (Fig.   3).     The 
bi -modal    character  is  less  marked   in  the  summer  months  when  there  is 
strong  lake  stratification.     The  longshore  currents  become  weaker, 
less  persistent,  and  are   interspersed  with  weak  or  rapidly 
fluctuating  current  flows  and  other  miscellaneous  events  (Fig.  4). 

The  following  current  regimes   have  been  identified  among  the  set  of 
miscellaneous  current  events. 

a)       Short  term  shore  parallel    and  antiparallel   current  episodes 
(~1  day). 


b)  Weak  onshore  and  offshore  current  episodes. 

c)  Low  speed  current  episodes  with  random  speed  and  direction 
fl uctuations. 

d)  Stagnation  episodes. 

These  current  flows  are  too  complex  and  irregular  to  be  analysed 
with  our  quasi-steady  state  model.     Due  to  their  low  speed  and  short 
duration  they  do  not  disperse  but  merely  pool    the  effluents  near  the 
source  (Murthy  and  Kenney  1974).     The  resulting  patch   is  not  only 
undesirable  from  the  viewpoint  of  effluent  dispersal    in  the  coastal 
zone,  but  under  unfavourable  conditions  the  pollutants  may  be 
transported  toward  and  trapped  near  the  shore  (Murthy  1972).     Such 
an   incident   is  all   the  more  likely  because  the  mass  exchange  of 
coastal    and  offshore  waters  preceding  the  development  of  a 
persistent   longshore  current  may  not  occur  between  successive  low  or 
random  current  episodes. 

4.      DIFFUSION  AND  THE  GAUSSIAN   PLUME 

The  spreading  of  a  dye  plume  in  a  steady  coastal   current  in  the  wake 
of  a  continous  source  exhibits  two  different  scales  of  growth  and 
diffusion   {Fig.   5).     The  first   is  the  lateral   growth  of  the  plume 
width   relative  to  the   instantaneous  centre  of  gravity.     This 
relative  diffusion  is  caused  by  turbulent  eddies  comparable  in  size 
to  the  plume  width   (Csanady  1963,  Murthy  and  Csanady  1971).     The 
second   is  the  meandering  of  the  plume  centre  of  gravity  caused  by 
turbulent  eddies  with  length  scales  larger  than  the  plume  width 
(Murthy  and   Kenney  1974,  Murthy  and  Miners   1980).     Combining  or 
superimposing  these  two  statistically  independent  processes  gives 
the  net  or  absolute  diffusion  of  the  plume  (Gifford   1959,  Bowden  et 
al.   1974). 

It  has  long  been  recognized  that  at  any  instant  the  concentration 
field  of  an  effluent  plume  is  likely  to  be  neither  Gaussian  nor 
symmetric.     This  is  hardly  surprising  considering  the  complexity  of 
the  coastal   flows,   influenced  as  they  are  by  wind  and  current 
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shear.     The  average  plume  concentration  field,  however,  can  be 
described  by  a  syrmietric  cross  plume  Gaussian  distribution.     This  is 
shown  by  figures  6  and  7,  which  diagram  48  successive  crossings  of  a 
dye  plume  at  a  constant  distance  from  the  source  (Murthy  and  Miners 
1978,   1980).     They  show  that  the  relative  cross  plume  concentration 
distribution   is  nearly  Gaussian  and  symmetric  (Fig.  6),  however,  the 
absolute  cross  plume  concentration  distribution  is  flatter  and  has  a 
sharper  cutoff  than  a  Gaussian  distribution  (Fig.  7). 

This  suggests  the  following  strategy  when  utilizing  time  series 
current  meter  data  to  model   the  average  concentration  field  of  an 
effluent  plume.     An  analytic  cross  plume  Gaussian  model    is  used  to 
calculate  the  concentration  field  relative  to  the  plume  centre  of 
gravity  for  each  different  speed  and  direction  event  of  a  shore 
parallel   episode.     These  realizations  are  then  weighted  and  averaged 
together  according  to  their  frequency  of  occurence  during  the 
current  episode.     The  final   numeric  -  analytic  estimate  of  the 
average  concentration  field  developed  by  the  selected   flow  regimes 
thus   includes,    in  an  approximate  manner,  the  combined  effects  of 
relative  diffusion,  current  speed  fluctuation,  and  plume  meandering. 

Having  discussed  the  formulations  of  the  advective  flow  field  for 
the  model,  attention  is  now  directed  to  the  diffusion  formulations. 
Three  different  diffusion  models  are  generally  considered  (Murthy 
1976).     There  is  Fickian  diffusion  where  the  eddy  diffusivity  is 
constant,   shear  diffusion  where  the  eddy  diffusivity  grows  linearly 
with  the  scale  of  the  diffusion  field,  and  inertial   subrange 
diffusion  where  the  eddy  diffusivity  has  a  4/3  power  dependence  on 
the  scale  of  the  diffusion  field.     Their  properties  are  sumnarized 
in  Table  1.     Now  the  scale  and  intensity  of  the  turbulent  eddies 
inside  the  coastal    zone  depend  not  only  on  the  local   wind  and 
current  flows  but  also  on  the  history  of  past  lakewide  wind  and 
current  events.     The  difficulty  is  to  decide  which  diffusion 
mechanism  or  combination  of  diffusion  mechanisms   is  operative  at  any 
time  at  a  given  site. 

The  following  conclusions  can  be  drawn  from  past  coastal    zone 
experiments.     The  vertical  diffusivity  is  at  least  one  to  two  orders 


of  magnitude  less  than  the  horizontal   diffusivity  and  can  be 
neglected   in  most  plume  spreading  calculations  (Csanady  1970,  Murthy 
1974).     The  difference  in  magnitude  is  attributed  to  stratification 
effects  and  the  difference  between  the  horizontal   and  vertical 
scales  of  motion  in  the  coastal    zone.     It  has  also  been  observed 
that   in  a  steady  coastal  current  the  width  of  an  effluent  plume 
increases  according  to  the  4/3  power  law  in  the  immediate  vicinity 
of  a  point  source  outfall    (Csanady  1970,  Murthy  and  Csanady  1971), 
however,  the  mechanisms  of  plume  growth  are  less  clear  at  distances 
far  from  the  source.     The  centreline  concentration  has  been  observed 
to  fall   off  as  x"1^2  (Csanady  1970,  Murthy  and  Csanady  1971),  and 
x       (Murthy   1972,   1974),  indicating  that  Fickian  and  vertical    and 
horizontal    shear  diffusion  processes  are  present.     Investigations  by 
Murthy  and  Miners   (1978),  Lam  and  Murthy  (1978)  suggest  that  the 
horizontal   eddy  diffusivity  of  the  coastal   zone  can  be  described  by 
the  empirical    relationship  shown  in  Figure  8,  however,  with  suitable 
coefficients  the  standard  diffusion  models  are  equally  applicable. 

In  summary,  the  ability  to  model    the  dispersion  of  effluents  in  the 
coastal    zone  is  directly  related  to  our  knowledge  and  understanding 
of  coastal    zone  diffusion  processes.     The  average  concentration 
field  of  an  effluent  plume  can  be  modelled  with  reasonable  accuracy 
if  site  specific  measurements  of  the  eddy  diffusivity  under 
different  flow  regimes  are  available.     Such  concentration  estimates 
are  less  accurate  if  the  strength  and  character  of  local   diffusion 
mechanisms  must  be  inferred  from  the  data  of  similar,  but  not 
necessarily  equivalent  coastal    sites,  or  from  diffusion  diagrams. 
(Okubo,   1971  and  1976) 

Kohl i    (1981)  reviewed  the  farfield  dilution  methods  with  steady 
current  and  presented  a   simple  and  quick  method  of  estimating  the 
two-dimensional   dilution  field.     This  method  may  be  used  for  quick 
and   rough  estimates  of  dilution  ratios  at  different  distances  from 
the  sources  under  steady  current  regime. 


5.  THE  GAUSSIAN  PLUME  MODEL 

The  analytic  equations  describing  the  diffusion  mechanisms  and 
average  concentration  field  of  an  effluent  plume  in  a  steady  current 
in  the  wake  of  a  line  source  diffuser  are  readily  available  in  the 
literature  (Fisher  et  al  1979,  Hinze  1959).  The  equations  and  their 
derivations  are  repeated  here  for  completeness  following  the 
analysis  of  Brooks  (I960).  The  advection  -  diffusion  equation  of  a 
non -conservative  tracer  with  an  exponential  decay  rate  in  a  field  of 
homogeneous  turbulence  is 


frbfrt^HH(K$H{*$-x< 


d) 


where  c(x,y,z)  is  the  tracer  concentration,  (u,v,w)  are  the 
velocities  in  the  (x.y  z)  directions,  (Kx>  K  ,  Kz)  are  the 
corresponding  eddy  diffusivities,  and  X  is  the  decay  lifetime  or 
dieoff  constant  of  the  tracer.  Consider  now  a  steady  and  continuous 
effluent  line  source  of  length  b  kept  perpendicular  to  a  uniform  and 
steady  shore  parallel  current  of  speed  U  as  shown  in  figure  9.  To 
allow  for  initial  jet  mixing  it  is  assumed  that  the  effluent  has 
been  diluted  to  concentration  c0  prior  to  release  from  the  ideal 
line  source.  Under  these  conditions  equation  (1)  reduces  to 
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(2) 


The  following  assumptions  further  simplify  the  analysis: 

1)  It  is  assumed  that  the  effluent  field  formed  by  the  line 
source  moves  downstream  in  the  x-di recti  on  at  the  same  rate 
as  the  prevailing  current  and  without  disturbing  the 
existing  flow  pattern  of  the  coastal  zone. 

2)  The  diffusion  in  the  flow  direction  is  negligible  in 
comparison  to  the  advection  (K x~o).  The  assumption  is 


justified  by  the  fact  that  the  concentration  gradients  in 
the  x-direction  are  very  small    compared  to  the  lateral 
direction  because  of  the  advection  in  the  x-direction. 

3)  Vertical   diffusion  is  negligible  compared  to  the  lateral 
horizontal   diffusion  (K    »KZ  ~  o).     As  described  in 

the  preceding  section  this  is  due  to  the  difference  between 
the  width  and  depth  scales  of  the  coastal   zone. 

4)  The  effluent   is  distributed  uniformly  over  the  available 
depth,  which  is  assumed  to  be  constant,  so  that  the  mixing 
can  be  described  by  a  two-dimensional  depth  integrated 
analysis.     The  available  depth  may  be  a  section  of  the 
water  column,  the  distance  from  the  thermocline  to  the 
bottom,  a  surface  layer,  or  the  total   height  of  the  water 
column,  depending  on  the  diffuser  characteristics  and 
climatological   conditions. 

5)  The  lateral  eddy  diffusivity  K    is  assumed  to  be  a 
function  of  L,  the  plume  width,  which  with  the  preceding 
assumptions  is  only  a  function  of  x,  the  distance  from  the 
source,   and  not  of  y.     The  assumption  implies  that 

Ky/K=g[L(x)]/g(b)=f(x) 

so  that  if  an  eddy  diffusivity  K  =  g(b)  is  associated  with 
the  length  b  of  the  diffuser  then 

K  =  g[l_(x)]=  f(x) 

where  L(x),  g(x),  and  f(x)  are  functional    relationships. 

6)  The  outfall    is  located  sufficiently  far  offshore  that  the 
spread  of  effluent  in  the  coastal    zone  is  not  restricted  by 
the  shore  boundary. 

With  the  above  assumptions  the  advection  -  diffusion 

equation  reduces  to 

&_K  dfc+x      0 

dx        ydy2  (3) 
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Equation   (3)  can  be   solved   in  the  half  plane     x^o 
the  boundary  conditions  of  x  -  o  being 


•I 


C  =  c0      for  |y|  <  b/2 

c  =  0         for  |yj  >  b/2  (4) 


The  dieoff  term  is  eliminated  by  the  variable  change 

cCe  (5) 


which  simplifies  equation   (3)   to  the  form 

K   £c  =  u2c 

df       dx   (6) 

where  C(x,y)  is  the  concentration  with  no  dieoff  and  x/u  is  the 
travel  time  to  any  distance  x.  Thus  the  diffusion  problem  can  be 
solved  without  dieoff  and  then  multiplied  by  the  decay  factor 
e"  x/u •  The  spatial  dependence  of  the  eddy  diffusivity 

is  eliminated  by  the  transformation  (assumption  5  above) 

dX  =  *<y  _  g(l)  _ 

dx   K   g(b)      (7) 

which  reduces  equation   (6)   to  the  form 

K^C  =  udc 
dy2  "     dx (8) 

The  solution  of  this  equation  with  standard  integral  techniques 
(Hildebrand  1965)  is  found  to  be 

c(x,y)  =  C(x,y)  e~Xx/u  =  C(X,y)  e~Xx/u  cga> 

b/2 


,-Xx/u 

/expJ-(y-tf)74KX/u{  dfl 

(9b) 
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\ar([   ty+b/2)1  -T   (y-b/2)l| 

)erl^Kx7rJ-erf[v^x7rr-J[ (9c) 
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Before  equation   (9)  can  be  utilized  the  functional    relationship 
between  X  and  x  must  be  specified   (equation  7).     This  is 
accomplished  as  follows.     Set 


L  s  2^/3  -a 


(11) 


where    a    is  the  standard  deviation  of  the  cross  plume  distribution 
function,  i.e. 

oo  oo 


^  =  ctb/y2c^>dv  =  cfe/yac(X^dy 


(12) 
-do  -co 

With  these  definitions  L=b  at  x=o  and  equation   (12)  can  be 
integrated  to  yield 

Qt"&f '[•{*?)&] 

Now  from  the  preceding  section  the  standard  diffusion  mechanisms 
relating  the  lateral   eddy 
of  the  effluent  plume  are 


(13) 


relating  the  lateral   eddy  diffusivity  K    to  the  scale  or  width  L 


1)  Fickian  diffusion  K  =  const  f(x)  -  1   (14a) 

2)  Shear  diffusion  K  a      I  f(x)  =  L/b    (14b) 

3)  Inertial   sub-range  K  a      L4/3  f(x)=(L/b)4/3  (14c) 

diffusion 


Substitute  equation   (13)  and  any  one  of  the  functional    relationships  of 
equation   (14)   into  equation   (7).      Integrating  and  rearranging  the  result 
then  gives 


4KX  =  b2jr      U4KA    6ln_ 
u  6  )[_     nV   u   /  b*J 


\    (15) 


where  the  diffusion  model    index  n  ■  1,2,3  for  Fickian,  lateral    shear, 
and   inertial    sub-range  diffusion  respectively.     The  results  are 
sunmarized   in  Table   1     and  Figure  9A. 
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The  most  useful   form  of  equation   (9c)  is  to  refer  all   concentrations  to 
the  effluent  concentration  c     bei 
Introducing  the  dilution  factors 


the  effluent  concentration  c     before  release  from  the  diffuser. 


TJ    =     §r     a  effluent  field  to  source  dilution   factor 

^?so=     c|     =  Jet  or  nearfield  to  source  dilution   factor 

allows  equation    (9c)  to  be  written 

71  -   ±n     c~Xx/u  j  crf  f  (y+b/2)1    crf|"   (y-b/2)"|  / 

'/(x.y)   -    2Vsoe  |ertLV4KX7^J     erfLV4KX7U-lj   (16) 

In  the  extreme  farfield,  the  line  source  approaches  a  point  source  and 
the  limit  of  equation   (16)   reduces  to 

71  -    Yt        p-*x/u    b      -y2/(4KX/u) 

7(x,y)  -  7so*e     V4nKX/u  '  e  (17) 

Equations  (15)-(17)  are  the  basic  equations  describing  the  average 
concentration  field  of  an  effluent  plume  in  a  steady  current  in  the  wake 
of  a  line  source  diffuser. 

6.  EXTENDED  GAUSSIAN  PLUME  MODEL 

The  Gaussian  plume  model  can  be  extended  to  include  the  effects  of 
current  speed  fluctuations  and  plume  meandering  in  the  final  estimate  of 
the  average  effluent  field  concentrations.  The  method  is  to  calculate 
the  equilibrium  or  steady  state  concentration  field  produced  by  all  the 
unique  speed  and  direction  events  of  a  shore  parallel  episode.  These 
realizations  are  then  weighted  and  averaged  together  according  to  the 
frequency  of  occurrence  of  each  generating  speed-direction  event.  Such 
a  procedure  presumes  that  the  average  distribution  of  concentrations  in 
a  quasi  steady  current  can  be  approximated  by  the  weighted  average  of 
successive  steady  state  solutions.  The  approximation  would  be 
reasonably  accurate  if  there  were  only  speed  and  no  direction  changes. 
This  is  because  the  above  weighted  average  does  not  take  into  account 
the  fact  that  the  fluctuations  in  the  flow  direction  of  a  meandering 
plume  tend  to  alternate  so  that  the  flow  remains  close  to  the  mean 
current  direction  (Fig. 5).  Consequently,  the  model  overestimates  the 
average  plume  width  and  underestimates  the  concentrations  near  the 
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centreline.     However,   it  has  been  observed  that  during  the  meandering  of 
a  plume  there   is  accelerated  mixing  due  to  shear  diffusion   (Murthy  1972, 
Murthy  and  Kenney  1974).     This  tends  to  compensate  the  limitations  in 
the  averaging  procedure  so  that  the  calculated  concentrations  are  likely 
to  be  better  estimates  of  the  average  field  concentrations  than 
originally  expected. 

The  averaging  procedure  is  applied  to  the  Gaussian  plume  model   as 
follows.     Consider  a  polar  coordinate  grid  (r,0)  centred  on  the  effluent 
source  with  the  principal    axis  parallel    to  the  shore   (Fig.    10).     Let 
there  be  a  current  u  passing  over  the  source  making  an  angle  0  with  the 
shoreline.     The  current  defines  the  x-y  coordinate  system  of  the 
Gaussian  plume  model   described  in  the  preceding  section. 

Koting  that 

x  =  rCos(0-<£)  y  =  rSin<0~0)    (18) 

and  the  restriction 


e-<t> 


<     n/2 


(19) 


the  dilution  factor  *7(x,y)  of  equation   (16)  can  be  written 

^(x.y)   =   *?(r,0,U.0)    (20) 

where  the  dependence  of  the  dilution   factor  on  the  speed  and  current 
direction  has  been  explicitly  included. 

Next  divide  the  time  series  data   (U,<£  )  of  a   shore  parallel    current 
episode  into  Q  segments  so  that  the  speed  and  direction   (U   .0   )  may 
be  considered  constant  during  the  segment  q,   say  hourly  averages. 
Define  the   speed -direct ion   frequency  distribution    fmn     of  the  current 
episode  by  sorting  the  speed-direction  segments   (U  ,  <f>   )  into  the 
M*N  bins  of  a  speed-direction  histogram  and  then  dividing  by  Q  so  that 

y^fmn  =  1 

™ (21) 
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This  establishes  the  statistics  of  climatology  of  the  current  episode. 
For  the  Great  Lakes  a  convenient   selection  of  histogram  intervals   is 
1  an/s  for  speed  and  5-45°  for  direction. 

Based  on  the  described   speed-direction  histogram  and  equations   (15)   - 
(21)  the  average  dilution  factor  at  any  (r,0)   is 

*?Cr.0)    =Ef'""7?Cr,0,Um,0„> 

«""  (22) 

An  alternate  method  of  calculating  the  average  dilution  factor  at  (r,0) 
is  to  first  calculate  the  average  sector  velocity   On    using  the 

equations 

In—    j    I  mn 

m  (23) 

On   =^U    Um/fm 

m  (24) 

and  then  substituting  the  result  in  equation  (15)  -  (21)  to  obtain 

n (25) 

For  well  behaved  current  episodes  the  difference  between  the  results  of 
equations  (22)  and  (25)  is  usually  not  significant. 

The  preceding  analysis  deals  with  a  single  current  episode  of  an 
identifiable  flow  regime,  say  a  shore  parallel,  medium  speed  flow.  If 
there  is  sufficient  data  then  it  is  possible  to  obtain  further 
cl  imatological  information  about  the  dispersion  and  mixing  properties  of 
the  selected  flow  regime  by  calculating  the  average  concentration  field 
of  each  current  episode  belonging  to  the  flow  regime  and  then  performing 
the  relevant  statistics  on  the  ensemble  of  average  concentrations. 

7.  MODEL  PARAMETERS 

The  input  parameters  of  the  extended  Gaussian  plume  model  are  those  of 
an  ideal  line  source  of  fixed  strength  lying  perpendicular  to  a  steady 
current  in  a  homogeneous  field  of  turbulence.  The  relationship  between 
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these  idealizations  and  their  physically  observable  equivalents  is  shown 
in  Figure  11.     It  shows  that  the  calculated  average  dilution  field  of  a 
line  source  diffuser  can  be  divided  into  a  nearfield  and  a  farfield 
solution.     The  nearfield  or  initial  dilution,  which  is  assumed  to  occur 
in  the  imnediate  vicinity  of  the  diffuser,  is  due  to  the  kinetic  and 
buoyant  potential   energy  of  the  diffuser  jet  discharge.     This  energy  is 
dissipated  in  the  turbulent  mixing  of  the  jet  with  the  surrounding  lake 
water.     The  degree  of  mixing  is  a  function  of  the  effluent  properties, 
water  depth,  diffuser  design  and  discharge  flow  characteristics.     After 
the   initial   jet  mixing  the  effluent  field  stabilizes  and  moves  with  the 
speed  of  the  ambient  field.     The  parameters  that  define  the  effluent 
field  or  nearfield  solution  at  this  juncture  are  the  ones  used  as   input 
parameters  to  the  farfield  or  Gaussian  plume  calculations.     They  are 
discussed  below  along  with  the  assumptions  and  simplifications   required 
to  estimate  their  values. 

7.1  The  Decay  Constant  or  Decay  Lifetime  (  t   =    1/X    ) 

The  inverse  of  the  decay  constant,  or  decay  lifetime,  is  the  time 
required  for  the  concentration  of  a  non-conservative  tracer,  obeying  an 
exponential   decay  law,  to  decrease  to  e~     of  the  initial   value.     The 
decay  constant  can   be  obtained  from  laboratory  measurements   if  the  decay 
rate  of  the  tracer  is  approximately  independent  of  the  ambient 
conditions  as   is  the  case  for  radionuclides.     The  experiments  are  more 
complex  for  those  chemical   and  biological   tracers  (e.g.  chlorine, 
coliform  bacteria,  etc.)  where  the  decay   rate  depends  on  ambient 
conditions  like  temperature,  mortality,   flocculation,  and 
sedimentation.      In  such  cases,    if  the  responsible  regulatory  agency  has 
not  set  the  decay  rates,   it  may  be  necessary  to  augment   laboratory 
results  with  on  site  measurements  at  different  seasons. 

The  importance  of  the  decay  lifetime  (t)  in  effluent  dilution  depends 
upon   its  magnitude  relative  to  the  time   (t)   required  to  achieve  dilution 
by  diffusion  alone.     If  r  »  t   then  the  first  exponential   term  in 
equations   (16)  and   (17)   is  nearly  unity  and  can  be  omitted,  i.e. 
e~AX/u   __   |  an(j  t^e  ef f i uent  can  be  treated  as  a  conservative 
substance  for  the  time  scale  of  the  calculations.     For  other  cases  the 
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complete  formulae  must  be  used.     However,   if  r«t   a  quick  estimate  of 
the  radial   extension  of  the  plume  for  a  desired  dilution       can  be 
obtained   from  the  formula. 

I       'so  e  ... (26) 

7.2  The  Diffusion  Coefficient   (K) 

The  selection  of  the  initial    value  and  growth  law  of  the  diffusion 
coefficient   in  the  coastal   zone  was  discussed  in  Section  4.     If  on  site 
measurements  are  available,  or   if  the  eddy  diffusivity  can  be  inferred 
from  the  measurements  at  similar  sites  with  similar  growth  laws,  then 
these  results  should  be  used.     Otherwise,  the  best  estimate  of  the 
diffusion  coefficient   in  the  coastal   zone  can  be  obtained  from  the 
empirical   growth  law  of  Figure  8  or  something  similar. 

7.3  The  Pi ff user  Length  (b) 

The  length  of  the  ideal    line  source  which  initiates  the  farfield 
calculations   is  the  width  of  the  effluent  field  after  the  initial 
dilution   (Fig.   11).     At  this  point  the  effluent  field  moves  with  the 
prevailing  current,   so  that  the  effective  line  source   is  perpendicular 
to  the  current,   regardless  of  the  original   diffuses,   -  current 
orientation.      If  the  current   is   roughly  perpendicular  to  the  diffuser, 
then  the  length  of  the  effective  line  source  is  about  the  same  as  that 
of  the  diffuser  projected  perpendicular  to  the  current  direction.     This 
approximation  clearly  fails  when  the  current  is  parallel   to  the  diffuser 
and  a  more   in  depth  analysis   is  then   required   (e.g.  Roberts  1977).      It 
is  also  assumed  that  the  distance  downstream  from  the  diffuser  to  the 
effective  line  source  is  small   and  can  be  ignored   in  comparison  to  the 
distances  of  the  farfield  calculations. 

7.4  The  Initial    Dilution  (  7}so  ) 

By  using  multiport  diffusers   it   is  possible  to   improve  the  dispersive 
efficiency  of  the  outfall   considerably. 
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The  initial  mixing  is  caused  by  turbulent  jet  diffusion,  a  flow 
phenomenon  that   is  governed  by  the  source  inputs  and  the  receiving  water 
characteristics.     A  clear  understanding  of  physical   processes  associated 
with  diffuser  outfall  discharges   is  needed  so  that  an  estimate  of 
initial   dilution  and  a  reliable  assessment  of  the  impact  of  wastewater 
discharges  on  the  environment  can  be  made. 

Hamdy  (1981)  reviewed  the  analyses  of  predictive  models  and  their 
applicability  to  different  types  of  submerged  diffuser  outfalls   in 
stagnant  receiving  waters. 

In  general,  for  the  purpose  of  initial   dilution  predictions  from 
submerged  sources,   a  distinction  has  to  be  made  between  "deep"  and 
"shallow"  receiving  water. 

Submerged  discharges  into  deep  receiving  waters  generally  exhibit  stable 
conditions  that   is,  upon  impingement  on  the  free  surface,  the  discharge 
spreads  laterally  in  the  form  of  a  stable  density  current.     These  stable 
conditions  create  a  stratified  near-field  flow  regime  which  allows  fresh 
anbient  waters  to  entrain  with  the  jet  flow  (Fig.   12).     The  vertical 
distance  below  the  water  surface  up  to  which  effective  entrainment 
occurs   is  dependent  on  the  thickness  of  the  surface  impingement  zone. 
The  diffusion  layer  thickness  is  the  depth  over  which  the  sewage  field 
is  uniformly  distributed  which  will   be  subjected  to  advection  and 
dispersion.     The  diffusion  layer  thickness  may  vary  from  8%  to  30%  of 
the  total   water  depth   (Abraham  1963).     In  shallow  receiving  waters, 
entrainment   into  the  jet   is  limited  and  the  jet  impingement  at  the  water 
surface  usually  results  in  a  complete  vertical   mixing  which  is  called 
unstable  near-field  (Fig.   13). 

Criteria  to  this  extent  have  been  developed  by  Jirka  and  Harleman 
(1973),  Almquist  and  Stolzenbach   (1976)  and  Lee,  Jirka  and  Harleman 
(1977). 

For  alternating  diffusers   (Fig.   14),  the  near-field  flow  regime  is 
unstable   if: 


H/D<  1.7  F  4/3   {i/D)'l/3 (27) 


18  - 


For  staged  or  unidirectional   diffusers   (Fig.   13),  the  near-field 
flow  is  unstable  if: 

H/D  <  0.46  F4/3   (      f/D)"1/3  (28) 

The  above  criteria  are  used  regardless  of  the  jet  angle      6     between 
the  jet  and  the  horizontal   plane. 

where 

H  =  water  depth  above  diffuser 

D  =  port  diameter 

F  ■  densi metric  Froude  Number  =  v/Vc^D 

n'      n*P 
9   =  9]D7 

V  =  jet  velocity 

Ap  *  the  density  difference  between  receiving  water 
and  eff  1  uent 

p»  =  the  receiving  water  density 

g  =  acceleration  due  to  gravity 

£  =  port  spacing  (between  two  consecutive  ports) 

7.4.1  Deep  Water  Models 

All  deep  water  models  are  based  on  an  integral  type  analysis  the 
principle  of  which  are  documented  by  Abraham  (1963),  Fan  (1967) 
Cederwall  (1971),  and  Fan  and  Brooks  (1969). 

Hamdy  (1981)  showed  that  these  models  may  be  limited  to  cases  of 
jets  with  no  jet  interference  and  with  low  discharge  densimetric 
Froude  numbers  (F«50).  The  above  limitation  is  valid  for  the 
shallow  coastal  zone  of  the  Great  Lakes. 
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Nomograms  shown  in  Figures  14  -  17  can  be  used  to  estimate  initial 
dilution  ratios  and  the  jet  width  for  horizontal  and  Inclined 
discharges. 

It  should  be  pointed  out  that  these  nomograms  predict  centreline 
ratios.     To  obtain  the  average  dilution,  multiply  the  centreline 
ratio  by  1.35  (Brooks  1973).     Also,  the  effective  water  depth  which 
is  about  70%  -  80%  of  the  total  water  depth  should  be  used  (Abraham 
1963). 

7.4.2  Shallow  Water  Models 

Models  for  different  diffuser  types  in  stagnant  shallow  receiving 
waters  predict  average  initial  dilution  in  the  nearfield  taking  into 
consideration  the  re-entrainment  effects. 

Estimates  of  average  initial  dilution  ratios  for  alternating  and 
staged  diff users  can  be  stated  as  follows  after  Jirka  and  Harleman 
(1973)   and  Almquist  and  Stolzenbach   (1976)   respectively. 


77^,(31  ternating)     =     .52   {  I  /  DF  )**    (H/D)    (29) 

77so(staged)  =  0.61  (th/D2ft  (30) 

Nomograms  shown  in  Figures  18  and  19,  can  be  used  to  easily 
determine  average  initial  dilution  ratios  for  alternating  and  stages 
diff users. 

7.4.3  Oiffusers  in  Ambient  Cross  Flow 

Roberts  (1977,  1979  and  1980)  has  reported  an  interesting  study  that 
reveals  important  three  dimensional   properties  of  line  diffusers  for 
buoyant  discharges.     Roberts'  experimental   study  and  dimensional 
analysis  for  currents  perpendicular  and  parallel   to  the  diffuser  are 
sunroarized   in  Figure  20. 
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The  estimates  are  valid  for  a  range  of  Froude  number  (F)  from  0.1  to 
100,  and  applicable  for  stable  near-field  flow  regimes.     For  sewage 
outfalls  with  negligable  density  differences,  the  near-field 
instabilities  and  associated  mixing  generated  by  the  momentum  flux 
tend  to  inhibit  any  formation  of  a  stably  stratified  zone  and  create 
a  uniformly  mixed  downstream. 

Thus,  in  the  case  of  unstable  near-field,  the  dilution  ratio  under 
cross  flow  will  amount  to  the  volume  of  flux  ratio  (  UHL/Q  ). 

8.      COMPUTER  SIMULATION 

A  computer  package  based  on  the  extended  Gaussian  plume  model   was 
developed.     It  has  a  section  for  calculating  the  coastal 
climatology,  the  average  concentration  field,  and  provisions  have 
been  made  for  users  to  read  in  the  results  of  section  7.4  or  to 
incorporate  their  own  initial   dilution  programs.     The  climatology 
section  sorts  hourly  average  wind  and  current  meter  data  into 
speed-direction  histograms,  prints,  plots,  and  stores  the  results 
for  later  retrieval   or  immediate  use.     The  speed-direction 
histograms  Are  used  in  the  diffusion  section  to  calculate,  print, 
plot,  and  store  the  average  concentration  fields  generated  during 
shore  parallel   current  episodes.     The  stored  concentration  fields  of 
several    shore  parallel   episodes  can  be  retrieved  and  combined  to 
obtain  the  average,  minimum,  and  maximum  concentrations  associated 
with  these  current  episodes. 

It  should  be  noted  that  what  is  desired  when  applying  the  above 
computer  model    is  not  the  average  concentration  field  calculated 
from  a  continuous  time  series  of  speed  and  direction  events  but 
rather  the  average  concentration  fields  associated  with  all   the 
individual    shore  parallel    and  antiparallel   current  episodes 
occurring  within  the  time  series.     The  concentration  fields  of  the 
individual   current  episodes  can  then,  under  suitable  conditions,  be 
averaged  or  combined  together.     The  necessity  for  this  approach  is 
illustrated  schematically  in  figures  21a-e. 
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Figure  21  shows  the  effluent  plume,  as  defined  by  a  dilution 
contour,  of  two  different   shore  parallel   current  episodes.     The 
current  speed   is  the   same  for  each  episode  but  their  directions  may 
differ.     In  Figure  21a  the  current  directions  are  the  same  so  that 
the  individual    plumes  and  their  average  are  the  same.     In  Figure  21b 
the  directions  of  the  two  plumes  differ  slightly,   as  might  occur 
when  a  plume  meanders.     The  average  plume  is  slightly  shorter  and 
broader  than  each   individual   plume.      In  Figure  21c  the  angle  between 
the  flow  directions  of  the  two  episodes  has  become  so  large  that 
their  average  concentration  field   is  only  half  as  long  but  somewhat 
wider  than  the  original    plumes.     In  Figures  21d  and  21e  it   is  clear 
that  the  outer  perimeter  of  the  two  plumes  (i.e.  the  concentration 
field  obtained  by  overlapping  the  two  plumes  and  retaining  the 
maximum  value)    is  a  better  estimate  of  the  combined  concentration 
field  for  our  purposes.     This  is  perhaps  the  main  difficulty  in 
using  the  extended  Gaussian  plume  model,   namely  deciding  how 
different  concentration  fields  should  be  usefully  combined. 

The  preceding   remarks  apply   in  a  somewhat   inverse  manner  to  the 
hourly  average  speed  and  direction  events  of  a  current   histogram. 
Suppose  a   histogram   is  obtained  from  two  current  measurements  having 
the  same  speed  but  different  directions.     If  the  angular  bin  width 
of  the  histogram  is  very  coarse  then  both  measurements  are  likely  to 
fall    in  the  same  direction  sector  leading  to  the  result  of  Figure  21a. 
Two  situations  can  occur  when  the  angular  grid   is  refined.      If  the 
angular  width  of  the  histogram  bins  or  the  angle  between  the  two 
currents  directions  is  small    compared  to  the  angular  plume  width, 
then  the  average  concentration  field  will   be  little  changed  as 
shown   in  Figure  21b.     However,    if  the  angular  difference   is 
significant,  then  the  average  concentration  field  will   be  reduced 
as  shown  in  Figure  21c.     This   indicates  that  the  original    sector 
angle  was  too  large.     It  should  be  decreased  until   the  sector  size 
no  longer  affects  the  calculation,   i.e.   until   the  diffusion  fields 
of  adjacent  current  sectors  overlap  as  shown  in  Figure  21b. 
Unfortunately,  this  condition  cannot  always  be  achieved  owing  to 
computer  memory  limitations. 
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The  four  January  current  episodes  of  Figure  2  are  now  analysed  as  an 
example.      It  will    be  shown  how  the  average  concentration  field  of  a 
single  current  episode  is  calculated  and  then  how  similar  episodes 
occurring  at  different  times  can  be  combined  to  give  an  estimate  of 
the  typical   dilution  field.     The  results  of  a  single  current  episode 
are  summarized   in  Tables  4-6  and   in  Figure  22.     Figure  22  shows  the 
plot  output  of  the  January  18-29  shore  antiparallel   episode  with  the 
100:1  dilution  contour  drawn  for  reference.     The  input  parameters 
for  this  simulation  were  the  diffusivity   (1000  cm2/s),  source 
dilution   (10:1),  diffuser  width   (200  m),  shoreline  angle   (250°), 
the  diffusivity  increasing  as  the  4/3  power  of  the  scale  of  the 
diffusion  field,  the  outer  radius   (10  km)  and  grid  sector  angle 
(5  ).     The  output  parameters  are  the  average  current  speed 
(11.8  cm/s)   and  the  diffusion  layer  thickness   (1.1  m),  which   is  less 
than  the  diffuser  depth  if  the  initial   dilution  estimates  are 
reasonable.     The  remaining  parameters  are  of  an   informative  nature 
and  not  used  in  the  calculations. 

The  current  advection  diagram  shows  that  the  current  carries  the 
effluent  to  the  edge  of  the  polar  grid,  as  it  must  if  steady  state 
calculations  are  to  apply  inside  this  region.     The  current  histogram 
or  rosette  shows  that  the  angular  resolution  was  fine  enough  to 
include  a  —15     plume  meander  in  the  calculations.     The 
speed-direction   histograms  of  Tables   2  and   3  show  that  this  episode 
consists  of  276  hourly  average  events,  that  there  is  a  plume  meander 
of   -15     about  the  mean  current  direction,  and  lists  a  maximum 
observed  speed  of  26  cm/s.     The  diffusion  field  of  Figure  22  is 
repeated   in   Table  4  which  lists  the  average  concentration  field  of 
this  episode  and  the  diffusion   field  overlap  between  adjacent 
current  histogram  sectors.     It  shows  that  the  diffusion  fields  of 
adjacent  current  sectors  do  not  overlap  after  the  first  3  km. 
Consequently,   the  sector  angles  should  be  decreased.     Since  this  is 
not  possible  due  to  computer  memory  limitations  the  calculation  will 
proceed  with  the  understanding  that  the  present  calculation  will 
overestimate  the  plume  length. 

The  calculations  were  carried  out   for  all    four  January  episodes. 
The  concentration  fields  of  the  two  shore  parallel   episodes  were 
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averaged  together,  as  were  those   of  the  two  shore  anti parallel 
episodes.     These   results  were  then  overlapped  and  combined  in  Figure 
23,  along  with  the  average  position  and  maximum  and  minimum  ranges 
of  the   100:1  dilution  contour. 

The  average  concentration  field  in  Figure  23  exhibits  a  bifurcate 
behaviour  both   in  the  shore  parallel   and  antiparallel   direction. 
This  may  be  characteristic  of  the  location  or  due  to  using  only  the 
four  January  episodes  instead  of  averaging  over  more  winter 
episodes.     Note  that  because  the  directions  were  different  it  was 
not   really  correct  to  average  either  the  shore  parallel   or  the 
antiparallel   episodes  together.     This  is  clearly  obvious  from  the 
maximum  limits  of  the   100:1  dilution  contour. 

9.     MODEL  APPLICATIONS 

Between  1970-72  the  Ontario  Ministry  of  the  Environment  gathered 
extensive  time   series  current  and  temperature  measurements   in  the 
vicinity  of  the  Lakeview  Water  Pollution  Control   Plant.     The  data 
were  used  to  assess  the  impact  on  adjacent  municipal  water  intakes 
of  different  Lakeview  WPCP  outfall    locations   (Kohli   and  Palmer 
1973).     The   purpose  of  this  section   is,   to  illustrate  how  the 
methodology  of  this  report  using  the   same  data  can  assist  in 
determining  dilution  contours  associated  with  this  outfall.     The 
procedure   is  to  calculate  the  average  dilution  field  of  a  particular 
effluent  from  several    realizations  of  a  given  current   regime.     The 
location  of  the   Lakeview  WPCP  and  current  meter  location  043  are 
shown   in  Figure   24.  The   shore  parallel   direction   is  taken  to  be  SW 
or  225  degrees  from  true  north. 

9.1  Data  Base  and  Procedure 

The  current  meter  operating  log  at  the  outfall    site   (043)   is 
presented   in  Table   5;  the  drogues  checked  reasonably  well   with  the 
current  meter  recordings.     The  data  are  the  pre-whitened  hourly 
average  current  speed  and  direction  measurements  obtained  from 
current  meter  location  043   (Kohli   and  Palmer  1973). 
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The  prewhitened  data  were  subjected  to  statistical    analysis  and  the 
summary  is  presented   in  Table  6.     The  currents  were  essentially 
bi -modal    in  character,  as  may  be  seen  from  the  histograms  of  Tables 
7  and  8.     The  time  series  data  was  divided   into  a  summer   (May  1  - 
Oct.   31)  and  a  winter  period   (Nov.    1  -  April    30).     Such  division 
showed  that  there  were  significant  differences   in  the  dimensions  of 
the  summer  and  winter  effluent  plumes.     Current  episodes  were 
selected   from  time  series  speed  and  direction  plots  and  stick  vector 
diagrams  for  each  of  the  two  seasons.     These  were  grouped   into  6 
current   regimes  according  to  three  speeds   (0-5,   5-15,  »15  cm/s) 
and  two  directions   (shore  parallel    and  antiparal lei ).     There  was 
some   subjecti veness   in  grouping  the  episodes.     For  example  the 
current  speed  of  an  episode  might  build  up  slowly,  then  shoot  up 
rapidly  to   12  cm/s  and   remain  there  a  while  before  returning  to 
zero,  the  average  speed  over  the  episode  being  4.9  cm/s.     The 
grouping  of  this  episode  then  depends  upon  the  relative  importance 
of  the  0-5  and  5-15  cm/s  contributions  of  the  episode.     With  minor 
exceptions,   the  average  current   speed  of  an  episode  fell   within   the 
range  of  the  current  regime  to  which  it  was  assigned.     (Tables  9  to   12) 
Episodes  of   less  than  a  day  were   rejected   in  order  to   insure  that 
there  was  a  well    developed  alongshore  current.     The  analysis  is  thus 
restricted  to  long  term  shore  parallel   transport.     The  effect  of 
short  term  or  non  shore  parallel   events  on  the  size  of  limited  use 
zones  must  be  obtained  by  other  methods. 

The  time-average  dilution  field  was  calculated  for  every  episode  of 
a  selected  current   regime.     These  realizations  were  then  combined  to 
produce  an  ensemble  minimum,  mean  and  maximum  dilution  field.     The 
resulting  dilution  fields  were  then  contoured.     For  a  given  current 
regime  this  procedure  not  only  provides  the  dimensions  of  the 
average  dilution  contour  but  also  an  estimate  of  its   range  of 
variation. 

The  results  of  the  average  dilution  contour  plots  are  summarized   in 
Table  13.     The  number  of  contour  plots  was  kept  at  a  tractable  level 
by  presenting  only  two  plots  for  each  current  regime,  i.e.  the  mean 
contours  and  the  minimum-mean-maximum  limits  of  a  selected  contour. 
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9.2   Input   Parameters 

The  input  parameters  used  for  the  computer  simulations  are  listed  in 
Table   14.     The  effective  length  of  the  diffuser,  i.e.  the  length 
projected  perpendicular  to  the  direction  of  current  flow,  was  set  at 
200m.     This  was  based  on  the  observation  that  the  average  angle  between 
the  normal   to  the  diffuser  axis  and  the  direction  of  current  flow  was 
about   5°  -  10°  for  well   developed  alongshore  currents  (Tables  7  and 
8),  accompanied  by  5°  -  15°  fluctuations  about  the  direction  of  mean 
current  flow  (Rosette  diagrams  of  Figures  1.01  -  1.20).     Thus,  the 
projected  diffuser  length  ranges  from  190-210  m  or  is  about  200  m  on 
average. 

The  selection  of  a  representative  diffusion  mechanism  for  computer 
simulations  is  somewhat  arbitrary  in  view  of  the  widely  differing 
diffusion  process  known  to  occur  in  the  coastal    zone  (Section  4).     An 
inertial-subrange  diffusion  model,  where  the  eddy  diffusivity  increases 
as  a  function  of  the  initial   plume  width  (Table  1),  was  adopted  because 
this  type  of  plume  dispersion  has  been  observed  to  occur  near  coastal 
zone  outfalls   (Csanady  1970,  Murthy  and  Csanady  1971).     The  predictions 
from  such  a  model  may  be  somewhat  optimistic  at  distances  far  from  the 
source  where  the  plume  width  becomes  comparable  to  the  offshore 
distance,  thereby  limiting  the  size  of  the  dispersing  eddies.     The 
diffusion  coefficient  associated  with  the  initial   plume  or  effective 
diffuser  width  was  therefore  chosen  very  conservatively.     The  diffusion 
coefficient  was  obtained  from  Figure  8,  divided  it  by  two,  and  rounded 
to  the  nearest   100  cm  /s.     Remembering  that  Figure  8  was  drawn  with 
the  convention  L  =  3a  ,  where    a    is  the  standard  deviation  of  the  cross 
plume  distribution  function,  while  our  model    uses  the  relationship  L  - 
2>/3-ct    (equation  11),  the  equation  providing  a  conservative  estimate  of 
the  eddy  diffusivity  near  the  line  source  becomes 

K     =     0. 01(^/2  •L)1"2 

»  0.01  (yJz/2* 20000) lmZ      where  L  =  200  m 

=  1200  an2/s 
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9.3     Estimate  of  Initial   Dilution 

The  following  outlines  physical   and  governing  parameters  for  the 
Lakeview  sewage  treatment  plant  outfall 


Discharge  Water       Port         Port  #  Jet  £/D  H/D 

Rate  (a)     Depth      Diam.       Spacing     ports     vel . 


m3/s  H(m)         D(m)  f(m)  m/s 


2.1  9.5  0.4  12  35         0.48  30  24 

The  discharge  densi metric  Froude  number  is  expressed  as: 

F  =  V/(g'D)   l/t 

Ap 
where  g'=  gpr 

For  an  effluent  temperature  of  22     and  the  receiving  water 
temperature  of  14     typical   temperatures  during  the  sunnier,     g'  = 
0.0146  and  F  =  6.3. 

As  the  port's  diameter  is  less  than  that  of  the  diffuser  outfall, 
the  effective  Froude  number  of  the  jet  is  usually  increased  by  a 
factor  of  about   1.8  (Effect  of  contraction  which  mainly  increases 
the  jet  velocity).     Thus,   the  effective  densimetric  Froude  number  is 
11.34. 

Since  the  diffuser  can  be  regarded  as  an  alternating  diffuser,  the 
near-field  flow  stability  is  examined  according  to  equation   (27) 
i.e.  H/D<  1.7  F4/3   (      f/D)   ~1/3. 

The  right  hand  of  the  above  equation  is  13.4  which  is  less  than  H/D 
indicating  a  stable  near-field  flow. 
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During  winter  conditions,  when  the  ambient  receiving  water  of  5°C 
and  an  effluent  temperature  of  15  c  may  prevail,  the  effective 
densimetric  Froude  number  is  about  20.     The  near-field  flow  regime 
in  this  case  may  be  considered  weakly  stable. 

For  stable  near-field  condition,  centreline  dilution  estimated  from 
nomogram  (Figure   14)   is   13:1   (an  effective  relative  water  depth  of 
0.8  x  24  m  20  is  considered)     The  average  dilution   is  13  x  1,35  = 
18. 

For  the  weakly  stable  condition,  it  is  appropriate  to  use  nomogram 
to  estimate  the  average  initial   dilution  ratio  which  will  amount  to 
16:1. 

The  above  estimates  are  for  stagnant  or  calm  conditions. 


Under  cross  flow  currents  the  initial   dilution  may  be  estimated 
according  to  Roberts'    (1979)  for  the  case  of  a  stable  nearfield. 

Based  on  a  current  speed  of  10  cm/s  perpendicular  to  the  diffuser 
axis,  an   initial   dilution   55:1  is  achieved.     As  the  current 
magnitude   increases,  the  initial   dilution  may  reach  100:1  or  more. 

For  the  weakly  stable  case,  the   initial   dilution  is  equivalent  to 
the  volume  flux   ratio,  i.e.         (    UHL/O   )   for  a  current  speed  of 
10  cm/s,  a  ratio  of  95:1  is  expected. 

The  calulations   in  this  appendix  use  an   initial   dilution  of  13:1  for 
all   current  flows.     This  conservative  estimate  was  obtained  by 
considering  the  stable  nearfield  centreline  dilution  under  stagnant 
conditions. 

9.4     Results  and  Discussions 

The  results  of  computer  simulations  based  on  the  extended  Gaussian 
plume  model    are  presented  in  Tables  7-13  and  the  contour  plots 
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1.01   -  1.20  as  described   in   Section  9.1.     Tables  7  and   8  summarize 
the   simmer  and  winter  coastal  current  climatology  near  the  Lakeview 
WPCP  in  the   form  of  a  current  speed-direction  histogram.     During 
both   seasons  the  currents  ape  either  mainly  shore  parallel    (233  deg. 
true)  or  shore  antiparallel    (38  deg.  true).     The  alongshore  current 
events  persisting  for  about  a  day  or  more  and   used  in  the  analysis 
are   listed   in  Tables   9-12  together  with  their  duration,  average 
velocity,   and  frequency  of  occurence,  the   latter  being   summarized   in 
Table   15.     The   frequency  Table   15  shows  that  the  present  analysis 
encompasses  42.2%  of  the   sumner  and  77.1%  of  the  winter  current 
measurements.     These   figures  would  be  even  higher  had  alongshore 
current  events   persisting  for  less  than  a  day  been   regularly 
included.     The  present  analysis   is,  therefore,  probably  valid  during 
most  of  the  winter  and  much  of  the   summer,   however,  there  will    be 
some  current  flows,  especially  during  the  summer,  that  will   disperse 
the  effluents  outside  the  limits  predicted  by  the  present  analysis. 

For  each  current  regime  the   reliability  of  the  dilution  contours 
depends  upon  the  number  of  episodes   used  to  calculate  them  and  their 
variability.     Thus  the  four  current  episodes  of  the  0-5  cm/s  speed 
regime  are   likely  to  produce  a   rougher  estimate  of  the  average 
dilution   field  than  the  eight  episodes  of  the   5-15  cm/s  speed 
regime.      In   such  cases   it   is  useful    to  combine  the  episodes  to  find 
the  maximum  and  minimum  limits  of  the  average  dilution   field    (see 
Table  3).     On  the  other  hand  the  very  high  current  episodes  are 
often   so  similar  to  each  other  that  the  concentration  calculated 
with  any  one  or  all   of  them  is  about  the   same   (e.g.   conservative 
effluent,  winter,  shore  parallel,   speed    >  15  cm/s). 

One  of  the   conclusions  that  emerges  from  an  examination  of  the 
results  is  that  there  is  a  significant  difference  between  the  summer 
and  winter  current  flows.     This  difference  manifests   itself  in 
several   ways.     The  summer  and  winter  speed-direction  histograms 
(Tables  7  and  8)   indicate  that  current  speeds  &re  smaller  and 
current  direction  fluctuations  larger  in  the   summer  than  in  the 
winter.     There   are   persistent  alongshore  currents  lasting   longer 
than  a  day  42.2%  of  the  time  in  the  summer  but  77.1%  of  the  time  in 
the  winter   (Table   15).     There  are  no  persistent  current  episodes 
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with  speeds    >  15  an/s  in  the  sunmer   (Tables  9  -  12).     Finally 
dilution   is  achieved  closer  to  the  source   in  sunnier  than   in  winter, 
possibly  because  the  current  meanders  more  in  the  sunmer. 

There  does  not  appear  to  be  a   significant  difference  between  the 
dilution  characteristics  of  shore  parallel   and  antiparallel    flows. 
In  the  summer  dilution  is  achieved  slightly  closer  to  the  source  by 
shore  parallel    than  by  shore  antiparallel    currents  belonging  to  the 
same  speed   range,   however,  the  converse  applies  to  the  winter  flows 
(Table   13).     For  a  given  speed   range  the  plume  size,  i.e.  the 
dimensions  of  an  average  concentration  contour,  depends  more  on  the 
season  than  the  current  direction. 

10.   FUTURE   DEVELOPMENTS 

This  background  report  has  described  the  coastal    zone,  some  of  the 
simpler  and  frequently  occurring  flow  regimes,   and  their  dispersion 
properties.     A  complete  description  of  the  coastal   currents,  their 
sources  and  driving  mechanisms,   is  not  yet  possible  but 
investigations  are  continuing  at   NWRI  towards  this  objective.     A 
computer  package  and  methodology  based  on     the  extended  Gaussian 
plume  model    has  been  presented.     The  model    package  using  a  limited 
number  of  time  series  current  meter  measurements  or  current 
climatology  data,   provides  an  approximate  but  rapid  estimate  of  the 
alongshore  effluent  transport.     Present  experience  suggests  two 
improvements  to  the  computer  model . 

1)  The  model    is  more  sensitive  to  the  current  direction  than 
the  current  speed.     The  accuracy  of  diffusion  field 
calculations  can  therefore  be  improved  by  resolving  the 
current  directions  to  the  nearest   1-2  degrees  and  the 
current   speeds  to  the  nearest  5  cm/s,  i.e.   the  reverse  of 
present  practice. 

2)  Short  duration  current  episodes,  where  the  plume  length 
is  limited  by  advection  rather  than  diffusion,   should  be 
included  in  the  analysis.     This  can  be  accomplished  by 
terminating  the  diffusion  field  calculation  of  these 
episodes  at  their  advective  transport  limit. 
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TABLE  1 

A  suimiary  of  the  diffusion  characteristics  and  lateral  spreading  of  an  effluent  plume  from  a  line  source  in  a 
steady  current. 
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TABLE  3 


CURRENT    SFEED-D1RECTIQN   HISTOGRAM         (TO   THE   NEAREST   TENTH   OF-    A   PERCENT) 

CONSERVED    (WINTER).      LAKE   ONTARIO   DATA.    PICKERING.       SHORE    ANT [PARALLEL    CURRENT. 

CURRENT   SECTORS         -      276  TOTAL 

IpFFn*      iNCRr^NT    "      ?'S  ^?lrS  SSS   B[N      '   H0LDS      «'»   CM'SEC   EI^ENTS.    ETC. 

SPEED         INCREMENT-       1.8CTVSEC.  SPEED   BIN   51    HOLDS      >  50  CIVSEC   EVENTS. 

SECTOR   CMRECTIONS  ARE   LABELLED   ANTICLOCKWISE   STARTING   FROM   THE   250   DEGREE   SHORELINE 

urMrcNIJo|NZ!ufi5^   ^Ei2HTED   ^6-RAGED   OVER   THE   NEAREST   NEIGHBOUR   SPEED  AND  DIRECTION  BINS 
HENCE    THERE   MAY   BE   SOME   ROUND   OFF    ERRORS    IN   THE    INTEGRAL   DISPLAY  OF    THE   H I STOGRAM  BELOW . 
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!■'  .0 
<*H.O 

n>.  0 

80.0  3  3      1 

85.0 

90.0  4  2     2 

95.0 

100.0  1 

105.0  10 

110.0  10 

115.0  12 

120.0  2 

125.0  26 

130.0  120 

135.0  193 

1-10.0  172 

145.0  109 

I50.O  '10 

155.0        re 
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I  »».«J 
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J29.e> 
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.'35.0 
2-10.0 
245.0 
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2-55.0 
260.0 
265.0 

**I?*2  ON  SHORE 

280.0 

205.0 

L-9D,0 

295.  a 

30(1.0 

30r-.0 
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515.0 

320.0 

3.25.0  4  )      g 

330.0 
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340.0 

345.0  6  2      3 

J50.0  1  1 
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TABLE  4 


COMPUTED  MEllN  CflNCIKTRfll  ION   I   II  LP 
(BffilO   ON  CIKRFN1    SPIFD  DIRECTION  HISTUGRAH1 


PARALLEL 

1.23 

KM 

OFF   SHQRf-: 

71.42 

KM 

ANT  (PARALLEL - 

8S.B2 

KM, 

ON     SHORE 

.39 

KM. 

""lCM  MWSEWES    (WINTtRl.      LAKE    ONTARIO  DATA.    PICKERING.      SHORE    AN1  .PARALLEL   CURRENT, 

F I  RSI    DAY      80'    118 

FINAL   DAY      80-    I/M  |,om,«  WWIMUH   COMPONENTS    OF 

.__    ._  "C   CURRENT   F'ECTOR   SUM 

DIFFUSION  CONSTANT  1000.00  CM**?' SEC 

SOURCE   DISCHARGE   RATE  2000   00  L'SEC 

DISCHARGE   VELOCITY  50.00  CM' SEC 

DENSITY   DIFFERENCE  [„S8  F>PT. 

SOURCE    DILUTION   RflTIO  10.0  (TO    1    0) 

MEAN   CURRENT   SPEED  u.8  CM^SEC 

DIFFUSION  LAYER  l.l  M    THICK. 

SIfFuIr     lInGT^  ,#2  fc                                                                                ORRENT  HISTOGRAM  DATA 

SHORElInE   SwU  MSI  clogr!                                                              CURRENT  "^CTORs""":"^"  TOTAL 

DIFFUSES    -SHORE   ANGLE  He  dEGREeI                                                              H££R        "SP*1"7    *      5.0      DEGREES 

DISTANCE      OFFSHORE  1.3  KM.                                                                         SPtED         '  CEMENT    -       1.0     CM^SEC. 

THE   DIFFUSIVITY    INCREASES   AS  THE   4/3   POWER  OF    THE  SCALE   OF   THE   DIFFUSION  FIELD 

CONCENTRATIONS   ARE  NORMALIZED   TO  SOURCE   CONCENTRATION  „TIM,Irn 

VALUES   PRESENTED  ARE    10..4   TIMES   ACTUAL   VHnis  CURREN? 

*?&£?  8gW  ieee-e    20ae-0     3eB0-e  40ae-0  5000.0     6000.0    7000.0     8e00.B    9000.0    10006.0  M.             mm 

™   B*         *      m      «:«  «:«  0-0      S-S      51      II      Si      s-2  P*^LEL    e-M 

10.0         240. 0                         .0                0.0                0.0  g;|  2-S                ?•?                8'0                |-«                0-0                0.0                                     0.00 

)'j.e      235.0               0.0           0.0           B.0  „'•  3  S'e           00           al           1*1           2-S           S-0                          0-a« 

20.0      230.0              0.0           0.a           0.0  J5  S/0           B-0           |.|           0-0          a.e           0.0                        e,00 

25.8        22S. 0                    0.0               0a               00  a  a  a   a               2'2               2'2               0-0               0-0               B-0                                  0-00 

38.0     22B.B                e          00          00  %■%  0-0          2-S          |-S          0-2         B-0          e-0                      0.08 

!S-0    ■?l->-0            -B        a  a        00  5:5  1/5        II        1*1        2-2       2-2        2-2                 0-00 

40.0        210.0                       .2                  .8                    B  "*2  °'S               *•■               g-g               B-0               0-0               0.0                                  0.00 

4b.8    265.0          2.0          s          1  :?  -S        0-2        0-0        B-B        a-2        e-e                 0-»« 

50.0   200.0        3.4      2.8      2.3  ,9  life      1.4      1  "I      «*f      ,  "f       *f             0-00 

55.0    195.0        2.3       .8        4  II  S        t      *•?      K       '-B       -9             0.00 

60.0    190.0         .5       .)        0  '§  *f       •"       ■'       -1       -1       .1             0.00 

65.0     185-e              .4                         *  t  a           'I           "I           "2           -2           -2                     0'00 

70.0        l80.0                    1.0                  .4                    3  'J  '1                  •%                  ••                  -8                  .0                  .0                                  0.00 

75.0         175.0                    2.0               1.0                  .7  ;s  ^                  -|                  ■>                  -J                  -1                  .1                                   0.00 


80.0         170.0                    3.0  1.6  1.1  .8  "e  S  'i  "I                  •?                  »l                                  0-00 

85.0         16b. 0                    3.S  1.9  U3  1  "I  -|  •?  -j                  'J                  .3                                  0.00 

3.5  »•*  ».a  •»  .'7  :i  "|  :3 

■"■"  1-7  1.1  .8  .7  6  S  3 

6.9  3.5  2.0  1.5  ill  ,:§  *f  ", 

16-7  8.8  6.2  4.9  4.0  3.4  3.~0  2  6 

20.5  11.9  8.5  u.t  5.6  4.7  4    1  37 

i^s    [30:5       Jf:f  ?:«  H  J-?  i*i  ?■?  f.j  ?.j       ;         i.z              -«•» 


<w.a  160.0  3.6  1.8  i.2  *l  '7  'I  "I  'i  -j  -8  0.00 

95.0  155.0  4.4  1.7  11  o  7  »|  -|  ■<  •"  -3  OFFSHORE           0.00 

1H0.0  150.0  8.9  3.5  2.0  u|  ill  .*!  ■!  "I  -4.  >i  0.00 

105.0  ,45.0  16.7  8.8  6.2  4.9  40  34  30  »"I  a  '-,  »"?  "-a" 

110.0  140.0  20.5  11.9  8.5  If  47  r?  S-f  H  ?"i  0-00 

US.0  135.0  18.3  8.7  5.8  4.3  3*1  »i  i*i  5"!  f5  ?.?.  0.00 


120.0  130.0  .39.2  7.9  a   ,  -,   ,  3-5  s-j  J"  '•'  i.»  J.f 

^0  IS::  ^:^  ,f:l  jg  i  #|  jg  »*./  ^  JM  ^^  ^;^ 

|:l  !!^:g  fltX  f£:J  J!§:^  !^:^  l  S  ffi  «•?  »  J^  M 

145.0  105.0  273.4  129.4  89  4  710  5n  B  S, '  3  15'q  ¥.-  »  *£-J  S6"8  112.15 

ill  IS  I;  II  K  I  I  I  II  1 1  P  B  « 

is  is  «  si  I  III  !  1  I  !l  sl 

ise.e  70.0  11.6  2.9  1.2  7  s  4  %  %  ■',  3-7,  0'00 

18-..0  66.6  1.5  .1  1"  -0  -0  "a  '»  'I  'I  '!      «NTIPARA           0.00 

lfl«-8  "0-0  •?  .0  0  ."  0  '0  *S  i  i  -2  B-eB 

195.0  55.0  2.1  .3  .o  2  "S  -0  -,B  -0  ■«  -8  0.00 

200.0  50.0  3.S  3.1  2.1  £4  Z'.t  I  "9  I-?  .*!  ,"0  ,-f  0-00 

205.0  45.0  2.4  .6  .3  a  p  I  *|  '*?  8*1  '*?  0-00 

I'IM.B  40.0  .3  .0  0  «  a  'a  'n  "i  •'  »'J  0-00 

1  1  II  li  N  H  h  I  1  P  8  § 

S»  351-0  0;a  S:S  lis  :  :  :  :  :  •■•  «  -0 

SSI  ^:l  5:5  g:^  S  ■  0  SI  :  It  I:S  "  -  —    5:31 

288.0  330.0  0.0  00  00  00  a'a  1*1  2'2  2-2  0-e  0-B  a-a0 

285.0  325.0  e.e  v.t  l\l  %■%  %-l  %■%  2-2  2-2  2-2  S-0  0-00 


•  ^a      0:0      1:1  1:1  1:1  1:5  %:%  %■%  %■%  ■■ 

0.0  0.0  0.0  0.0  0.0  0.0  h  e 

0.0  0.0  0.0  0.0  0.0  0.0  a'0 
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i!'..0  jr,.H  .8  1  a  *  a  "I  ■,  -2  -°  ■!  B-BB 

340.0  270.0  3:9  1  1  I  -2  -2  ■?  -B  -B  -°  a.«e 

345.0  26,;. 0  6.4  4.9  4.-0  J.4  2.'9  ,-|  *"j  ,  *J  ,  ' I  . -j  *>■*» 

350.0  260.0  4.6  1.9  1.1  .8  6  I  3*I  **!  l*f  Ul  0-00 

J55-«  «5-0  I J  -2  -1  0  ^  '.I  ".0  ".0  .8  ^ 
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RELATH.E  CENTRELINE  DECAY  AND  DIFFUSION  AT  MEAN  CURRENT  SPEED  INTO  AN  ADJACENT  CURRENT  OR  DIFFUSION  FIELD  SECTOR  | 100*  SOURCE. 

!:l       81     M     II     4:]     «:S     4     4     --     «•?     -d  ^L,^croR 

61-3  13-8  3-6  '-J  •<  -3  -2  -1  -1  .1      DIFFUSION   SECTOR 
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TABLE   5 
Location  043  Operating  Log 
Location:  79°  31'   54MW;  43°  34'   08"N 

Depth  of  Water:  31.0  feet  (9.5  m.) 

Depth  of  Current  Meter  9.8  feet  (3.0  m)  from  bottom 

Current  Meter  Details:  Serial   #  144  and  272,  Plessey  M021 


Date      Time  (Hours) 

In    Out  Service 


May  21  1409  Meter  installed 

Jun  8  1020   1005     Meter  checked  and  found  O.K. 

Jun  12  Drogue  trackings 

June  25  Drogue  trackings 

Jul  2  Drogue  trackings 

Jul  21  0945   0955     Meter  working  O.K. 

Jul  30  1210   1200     Meter  O.K.,  drogue  trackings 

Aug  13  1235   1220     Meter  changed,  drogue  trackings 

Sep  30  1240   1155     Changed  bearings,  "0"  ring,  battery, 

dessicator  and  tape.  Drogue  trackings 

(1971) 

DeTT"  1440  Meter  installed 

Dec   15  Surveyed  spar 


(1972) 
Fe 


Jeb  9  Drogue  trackings 

Mar  27  1155  Meter  removed 


-40- 
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TABLE   6 
Sumnary  of  Current  Meter  Results,   in  the  Vicinity  of  Lakeview  WPCP  outfall 

Jun   70     Jul    70     Aug   70     Sep   70     Oct   70     Nov   70     Dec   71     Jan    72     Feb   72     Mar   72 


Result  and  direction   100  109  176  60  189  68  62  56  57  223 

(0°       North) 

Resultant  speed 

(on/s)  1.38         0.34         1.24         1.06         0.83         2.31         1.81         5.43         0.44         1.24 

Average  speed 

(cm/s)  3.49         2.30         5.39         4.36         4.26         6.06         6.91         7.51         4.55         3.62 

Maximum  speed 

(cm/s)  20     14     20     20     17      27      34     29     24     29 

Persistence  factor   0.39    0.15    0.23    0.24    0.19    0.38    0.26    0.72    0.10    0.34 

Percentage  of  time 
going  in  direction 
of  resultant         4      2      4      12      1      13     23     41      18      19 

Percentage  of  time 

going  towards  R.  L. 

Clark  Intake         5     10      6     12     13     25     23     41     18     12 

Total  number  of 

readings  4312    4007    4461    4311    4464    4317    1584    2232    2088    1986 

Interval  of 

readings  (min)       10     10     10     10     10     10     20     20     20     20 

« 

Percentage  of 

negligible*  speed 

{%  of  recording 

period)  14  25  8  15  10  8  7  9  32  38 


TABLE  7 


CURRENT    SPEED- DIRECTION   HISTOGRAM         (TO    THE    Nf.rtRFPl    TENTH  OF    A  PERCENT) 

CONSERVED    i SUMMER).      ALL    MOE   SUMMER   LAKEVIEW  DAIA.      MAY      1 .  19 TO    -    OCT      31,l'=>?0. 


i.URRENT   VECTORS         -    3822    TOTAL 
SECTOR      INCREMENT    -      5.0  DEGREES 
SPEED         INCREMENT    -      1.0  CM'SEC. 


SPEED  BIN      1   HOLDS     0-1   CM'SEC  EVENTS.    ETC. 
SPEED  BIN  51   HOLDS      >50  CM--SEC  EULNTS. 


SECTOR  DIRECTIONS  ARE   LABELLED  ANTICLOCKWISE   STARTING  FROM    THE  22S  DEGREE  SHORELINE 

CURRENT   EVENTS  ARE  WEIGHTED  AVERAGED  OVER   THE   NEAREST  NEIGHBOUR   SPEED  AND  DIRECTION  BINS. 
HENCE    THERE   MAY   BE    SOME   ROUND  OFF   ERRORS    IN    THE    INTEGRAL   DISPLAY  OF    THE   HISTOGRAM  BELOW. 


SPEED 
SECTOR    i 

0.0 

5.0 

10.0 

15.0 
20.0 
25.0 
30.0 
35.0 
40.0 
15.0 
-.0.0 
55.0 

60.0 

65,0 

■  0.0 

75. e 

80.0 
85.0 

90.0 

95.0 

100.0 

105.0 

110.0 

115.0 
120.0 
125.0 
130.0 
135.0 
140.0 
145.0 
150.0 
155.0 

Ir.l.il 

165.0 

1  1*0.0 

175.0 
180.0 
185.0 
190.0 
195.0 
200.0 
205.0 
210.0 
215.0 
220.0 
295.0 
230.0 
235.0 
.'40.0 
24'...  0 
.150.0 
255.0 
260.0 
265.0 
2?0.0 
275.0 
280.0 
285. 0 
,.-90.0 
295.0 
300.0 
305.0 
310.0 
315.0 
320.0 
325.0 
330.0 
335.0 
340.0 
345.0 
350.0 
355.0 


BIN    • 
SUM 

28 
18 
17 
14 

14 
13 
II 

7 

9 
10 

7 

7 

7 

t 

7 

■< 

i 

6 

9 

B 

B 
II 

7 

9 

B 
1  I 
10 
13 
12 
12 
12 
12 
16 
16 
17 
30 
40 

51 
54 

51 
37 
29 
19 
13 
7 
2 
3 
2 
3 
5 
4 
4 
4 
4 
6 
5 
6 
B 
7 
9 
T 
9 
8 

!  I 
I  1 
11 
14 
19 
23 
27 
55 

53 


9      10      II      12      13      14      15      16      17      18      19     20     21 


■1 
5 
3 
4 
4 
5 
4 
3 
3 
5 
4 
4 
p 

i 

4 
5 
I 

2 

■! 

4 
4 

6 
! 
4 

4 
5 

4 
7 

7 
4 
■I 
4 
■I 
4 
•1 
9 
9 
9 
9 
9 
5 
i 

4 
4 
I 
1 

2 


2 
1 

2 
1 

1 
3 
2 
2 
2 
3 
2 
<■ 

2 
4 
3 
4 
4 
5 
3 


2 
1 
2 

2 

a 

2 

i 
i 

2 
2 
1 
: 
2 

I 

1 

2 
2 

a 

i 

2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
I 

2 
3 
3 
1 
4 
A 
3 

5 
4 
4 
4 


2 
1 

2 
I 
2 

2 
1 
1 
1 

I 

; 
i 
l 
i 
i 
i 
l 
l 

i 
i 
i 
i 
l 
I 
l 
t 
I 
i 
i 

2 

2 
3 
2 
2 
2 
3 
3 
5 
6 
5 
4 
3 
I 
I 


1 
I 
1 

1 
1 
1 
1 
I 

2 
1 

2 
2 
2 
4 
4 
7 
7 
5 
2 
1 
I 


1 

2 
1 
1 
1 
1 
1 
1 

2 
2 

2 
3 

2 
3 

2 


1 
1 
1 

l 
I 
1 
1 

2 
1 
2 
1 
2 
3 
5 
6 
5 
4 
3 
1 


1 
1 
1 

I 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
4 
3 


2 
A 
5 
4 
4 
1 
i 


1 
l 
l 
1 
2 
A 
5 
4 
3 
1 
1 


PARALLEL 


OFFSHORE 


? 

1 

1 
3 

A 

4 

4 
1 


ANT  I PARA 


ON  SHORE 


TOTAL 


1O00     270    123    110     94     84     68     57     53     3?     28     22     20      12 
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TABLE  8 


CURREN1    STEED- DIRECT  I  ON  HISTOGRAM        HO   TMfi   Nt.AKfSl    ILNTM  hi     fl  rtMlNl  i 

CONSERVED    (WINTER).      ALL   MOE   WINTER  LAkO'lEW  DAIA.      NOl'.     1,10ft)        ffl*RU     (1*19/ 1. 


t  URRENT   SECTORS         -    1958    TOTAL 
SECTOR       INCREMENT    -      5.0   DEGREES 
SPEED         INCREMENT    -       1.0   CM'SEC. 


SFEED   BIN      1    HOLDS      0    I    CM   SEC   El'tNIS,    tTt 
SPEED   BIN  51   HOLDS      >50   CM- SEC   EVENTS. 


SECTOR   DIRECTIONS   ARE    LABELLED   ANTICLOCKWISE   STARTING   FROM    THE   225   DEGREE   SHORELINE 

CURRENT   EVENTS   ARE   WEIGHTED  AVERAGED   OU-ER    THE   NEAREST  NEIGHBOUR   SPEED   AND  DIRECTION  BINS. 
HENCE   THERE   MAY   BE    SOME   ROUND  OFF   ERRORS    IN   THE    INTEGRAL   DISPLAY   OF   THE   HISTOGRAM  BELOW. 


SPEED 

SF.f  TOR 

0.B 
S.U 

10.0 
15.0 
20.0 
2b.  0 
30.0 
35.0 

40.0 

45.0 
50.0 
&S..8 

l«.0 
<V>.0 

rn .  0 
rs.a 

U0.0 
!.<5.0 
90.0 
9b.  0 

100.0 
100.0 
I  10.0 
115.0 

120.0 
125.0 
1  50.0 
135.0 

1-10.0 

1-15.0 

I '.,0.0 
155.  0 
100.0 
K.5.0 
170.0 

175.0 

130.0 

185.0 
190.0 
195.0 

200.0 

205.0 
210.0 
215.0 
.  20.0 
.V..0 
..'30.0 
235.0 
2-10,0 
2-15,0 
..•50.0 
,-55.0 
,.'60.0 
265.0 
270.0 
275.0 
280.0 
285.0 

290.0 

295.0 
300.0 
305.0 
510.0 
315.0 
V0.0 
325.0 
HI1.9 
'■35.0 
1-10,0 
'-15.  H 
'.0.0 
..55.0 


BIN    »    1 

SUM      -- 


39 
21 
12 
7 
6 
6 
3 
3 
6 


3 

3 

4 
7 

c, 

6 

5 

7 

E, 

(" 

12 

15 

22 

31 

4e 

57 
69 
53 
47 
28 
MS 
6 
5 


1 
2 
2 
2 
2 
1 
2 
2 
2 
2 
2 
1 

4 
■1 

A 

3 
7 
5 
6 

11 

1  15 

1  M'J 


2      3      4 


1 

1 
2 
2 
1 

3 
6 
4 
4 
3 
3 
1 
1 

6 
5 
6 
3 
6 
6 
A 
2 
2 
3 

1 

1 
1 


1       1 


1  1 

1  2 

2  2 
4  3 


1 

2 
3 
6 

10 

2 

S 
2 

1 


8     9    10    11    12    13    14    15    16   17    18    19  20  21   22  23  24  25  26  27   20   29   30 


1      1 


PARALLEL 


OFFSHORE 


ANT  I  PARA 


1        1       1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1    1 


1     1 


ON    SHORE 


10    1  1 
0    10 


2      3      111 


IOTmL       1O00    167   89   98  94    79    73  63  56  47   41    28  23   25    19    13    13    13     8      9    10      5      3      6      3      3      4 
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TABLE   9 

Summer  Current  Episodes  (Shore  Parallel) 
(Lakeview  U.P.C.P.) 


Speed  0-5  cm/s       (frequency  8.5%) 
Episode  Time  and  Date 

18:00     8/06/70  -  23:00  11/06/70 

12:00     4/09/70  -  23:00  7/09/70 

6:00     1/10/70  -  17:00  5/10/70 

12:00  11/10/70  -  17:00  13/10/70 

Average  episode 


Duration  Average  Velocity 

(hours)  (cm/s  &  deg.  true) 

78  1.9  at  245° 

84  4.5  at  226° 

108  2.5  at  223° 

54  4.2  at  221° 

40.5  3.1  at  2270 


Speed  5-15  cm/s       (frequency  11.0%) 


Epi  sode  Ti  me  and 

Date 

Duration 

Average  Velocity 

(hours) 
96 

(cm/ 

s  &  deg.   true) 

0:00  22/07/70  - 

23:00  25/08/70 

5.1  at  240O 

0:00  27/08/70  - 

11:00  29/08/70 

60 

7.1  at  230O 

12:00     1/09/70  - 

17:00     3/09/70 

54 

5.8  at  235° 

0:00  11/09/70  - 

17:00  12/09/70 

42 

8.6  at  2310 

0:00  17/09/70  - 

5:00  18/09/70 

30 

8.1  at  236° 

6:00  28/09/70  - 

5:00  30/09/70 

48 

4.3  at  230O 

0:00  17/10/70  - 

11:00  19/10/70 

60 

6.2  at  2320 

12:00  23/10/70  - 

17:00  24/10/70 

30 

6.4  at  230O 

Average  episode 

52.5 

6.2  at  233° 

Speed         15  cm/s 

(frequency  0%) 

No  episodes  observed 
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TABLE   10 


Sunnier  Current  Episodes  (Shore  Antiparallel ) 
(lakeview  W.P.C.P.) 


Speed  0-5  cm/s       (frequency  8.9%) 


Episode  Time  and 

Date 

Duration 
(hours) 

0:00  28/05/70  - 

23:00     4/06/70 

186 

0:00  18/08/70  - 

17:00  19/08/70 

42 

0:00     9/10/70  - 

11:00  11/10/70 

60 

12:00  29/10/70  - 

17:00  31/10/70 

54 

Average  episode 

85.5 

Speed   5-15  cm/s 

(frequency   13.8%] 

Episode  Time  and 

Date 

Duration 
(hours) 

12:00  24/05/70  - 

17:00  25/05/70 

30 

12:00  14/06/70  - 

11:00  17/06/70 

72 

0:00     8/08/70  - 

23:00  10/08/70 

72 

0:00  22/08/70  - 

23:00  23/08/70 

48 

6:00     8/09/70  - 

17:00  10/09/70 

60 

6:00  14/09/70  - 

5:00  16/09/70 

48 

18:00  19/10/70  - 

11:00  23/10/70 

90 

0:00  25/10/70  - 

11:00  29/10/70 

108 

Average  episode 

66.0 

Speed         15  cm/s 

(frequency  0%) 

Average  Velocity 
(cm/s  &  deg.   true) 

1.9  at  460 

4.0  at  380 

3.1  at  44° 
3.7  at  24° 

2.6  at  390 


Average  Velocity 
(cm/s  &  deg.   true) 

7.5  at  40° 

7.4  at  490 
8.0  at  430 
4.7  at  430 

4.5  at  420 
8.7  at  410 
4.9  at  430 
6.0  at  31° 

6.3  at  410 


No  episodes  observed 
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TABLE  11 

Winter  Current  Episodes  (Shore  Parallel) 

ILakeview  W.P.U.P.J 


Speed  0-5  cm/s       (frequency  10.0%) 
Episode  Time  and  Date 


12:00  29/11/70 

18:00  18/02/71 

18:00     8/03/71 

6:00  31/03/71 

Average  episode 


5:00  1/12/70 
23:00  20/02/71 
23:00    9/03/71 

5:00     3/04/71 


Duration 
(hours) 

41 
54 
30 
72 

49.3 


Average  Velocity 
(cm/s  &  deg.  true) 

3.2  at  2240 

3.3  at  238° 

5.2  at  240° 

4.3  at  233° 

3.9  at  234° 


Speed   5-15  cm/s       (frequency  25.7%) 


Epi  sode  Time  and 

Dat 

:e 

Duration 

Average  Velocity 

(hours) 
30 

(cm/s  &  deg.  true) 
10.3  at  2350 

06:00     6/11/70  - 

11:00     7/11/70 

12:00  16/11/70  - 

23 

:00  21/11/70 

138 

7.4  at  2340 

00:00  25/11/70  - 

11 

00  27/11/70 

60 

10.4  at  2340 

00:00     2/12/70  - 

17 

00     3/12/70 

42 

5.7  at  2340 

06:00     5/12/70  - 

11 

00     6/12/70 

30 

10.5  at  2410 

12:00     8/12/70  - 

11 

00  10/12/70 

48 

5.3  at  232° 

18:00  25/02/71   - 

17 

00  28/02/71 

71 

6.4  at  233° 

18:00     4/04/71  - 

05 

00     8/04/71 

84 

7.9  at  230O 

Average  episode 

63.0 

7.7  at  2340 

Speed       15  cm/s  (frequency  7.0%) 

Episode  Time  and  Date 


06:00  22/11/70  -  03:00  24/11/70 
06:00  2/03/71  -  23:00  4/03/71 
14:00  11/03/71  -  14:00  12/03/71 

Average  episode 


Duration  Average  Velocity 

(hours)  (cm/s  &  deg.   true) 

46  15.4  at  235° 

66  14.8  at  2350 

25  16.7  at  235° 

45.7  15.3  at  235° 
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TABLE   12 

Winter  Current  Episodes   (Shore  Antiparallel  ) 

(Lakeview  W.P.C.P.) 


Speed  0-5  cm/s       (frequency  12.3%) 
Episode  Time  and  Date 


6:00  11/11/70  -  23:00  15/11/70 
12:00  6/12/70  -  5:00  8/12/70 
12:00  29/03/71  -  5:00  31/03/71 
18:00     9/04/71  -  11:00  11/04/71 

Average  episode 


Duration 
(hours) 

Average  Velocity 
(cm/s  &  deg.  true) 

114 
42 
42 

42 

3.2  at  43° 
3.1  at  51° 
2.8  at  350 
1.5  at  61° 

60.0 

2.8  at  45° 

Speed       5-15  cm/s 
Episode  Time  and  Date 


(frequency  17.2%) 


12:00     2/11/70 
18:00     8/11/70 

0:00  28/11/70  ■ 
18:00     3/12/70 

0:00     1/03/71  • 

0:00  6/03/71  • 
12:00  26/03/71 

6:00  3/04/71  ■ 

Average  episode 


17:00     3/11/70 

5:00  10/11/70 

11:00  29/11/70 

5:00     5/12/70 

5:00     2/03/71 

11:00     8/03/71 

11:00  29/03/71 

17:00     4/04/71 


Duration 
(hours) 

30 
36 
36 
36 
30 
60 
72 
36 

42.0 


Average  Velocity 
(cm/s  &  deg.  true) 

8.5  at  30 
6.5  at  38° 
6.7  at  410 

7.1  at  320 

7.2  at  410 

5.3  at  55° 

5.4  at  42° 
6.3  at  390 

6.3  at  41° 


Speed       15  cm/s 
Episode  Time  and  Date 


0:00  21/02/71  ■ 
10:00  24/02/71  ■ 
15:00  10/03/71   ■ 

Average  episode 


(frequency  4.9%) 


17:00  22/02/71 
17:00  25/02/71 
12:00  11/03/71 


Duration 
(hours) 

42 
31 
22 

31.7 


Average  Velocity 
(cm/s  &  deg.   true) 


10.4  at 
16.6  at 

15.5  at 


30O 
350 
290 


13.6  at  320 
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TABLE    13 
Plume  Dimensions  at  the  Lakeview  W.P.C.P 

The  plume  dimensions  are  the  length  and  width  of  a  specified  dilution  contour  averaged  over  several  episodes.     The  maximum  and 
minimum  ranges  of  the  contour  are  also  given.     The  diffuser,  effluent,  and  current  episode  descriptors  are  listed  in  Table   15, 
the  most  useful  ones  being  repeated  below.     They  are  the  specified  dilution  (s),  the   frequency  of  occurrence,  and  the  diffusion 
time   (t)  required  to  achieve  the  specified  dilution. 


S     -     100  :   1 


SUKMER 


00 


Speed 
(cm/s) 

0-5 
5-15 
15 


Frequency 

m 

8.5 

11.0 


Shore 
Diffusion  Length  (km) 

Time  (min)     min      mean      max 


805 
670 


1.2 

1.0 


1.5 
2.5 


Parallel  Flow 

Width  (km) 

min       mean       max 


2.4 
6.0 


0.4 
0.3 


0.6 
0.5 


0.9 
1.0 


Frequency 

55 

8.9 
13.8 


Shore^  Antiparallel   Flow 
Diffusion  Length  (km)  Width  (km) 

Time(min)        min      mean      max        min      mean      max 


1155 
900 


1.3 
1.3 


1.8 
3.4 


4.0 
7.4 


0.3 
0.3 


0.6 
0.5 


0.9 
1.3 


WINTER 
Shore  Parallel   Flow 
Speed     Frequency    Diffusion  Length(kmT  Width  (km)  Frequency 

(cm/s)  (%)  Time  (min)     min      mean      max     min      mean      max  (%) 


0-5 
5-15 
15 


10.0 

25.7 

7.0 


1370 
1340 
1765 


1.0 
1.1 
1.0 


12.3 

17.2 

4.9 


Diffusion 
Time(min) 

1665 
1215 
1190 


Shore  Anti parallel  Flow 
Length   (kmT  Width  (km) 


mm       mean 


max 


min 


mean 


1.3 
1.3 

4.0 


2.8 
4.6 
9.7 


3.7 

7.1 

14.6 


0.3 
0.3 
0.5 


0.8 
0.5 
0.5 


max 


1.5 
1.5 
0.9 


TABLE   14 
Input  Parameters  for  Effluent  Plume  Simulations 
Initial    Dilution       -     13:1   {as  per  section  9.3) 
Diffuser  length         -     200  m 

Diffuser  orientation  -     135  deg.  true  (axis  direction) 

Shore  parallel  direction  -     225  deg,   true 

Eddy  diffusivity       -     1200  cm2/s. 


The  computer  simulations  are  based  on  the  inertial    sub-range  diffusion 
model,  i.e.   the  eddy  diffusivity  increases  as  the  4/3  power  of  the  plume 
width.     The  diffusion  constant  K  =  1200  cm^/sec  corresponds  to  an   initial 
plume  width  of  200  m. 


TABLE  B5 


Source  Description 


Discharge  rate  2100  1/sec  (typical) 

Discharge  velocity  48  cm/s                " 

Density  anomaly  1.5  ppt.               " 

Initial  dilution  13:1 


Diffuser  Depth  9.5  m 

Diffuser  Length  210  m 

Diffuser  Orientation  135  deg.  true  (axis  direction) 


Shoreline  Orientation         225  deg.  true 

Distance  Offshore  1.3  km 

Diffuser  Location  As  per  Figure  Al 
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TABLE   15 
Alongshore  Current   Frequency 

Tab 

e 

Lakevtew  W.P.C. 

Sumner  (May  1  -  Oct  31) 

(based  on  3822  hourly  averaged  events) 

Speed                                Parallel 
(cm/s)                                   (X) 

P. 

Ant  i  pa  rail  el 
(X) 

0-5                                       8.5 
5-15                                   11.0 
15 

8.9 
13.8 

Sum  19.5  22.7 

Percentage  of  alongshore  current  episodes  ■  42.2 
Percentage  of  miscellaneous  current  events  =  57.8 


Winter  (Nov  1  -  April  30) 

(based  on  1958  hourly  averaged  events) 

Speed              Parallel  Antiparallel 

(cm/s)              (X)  (X) 

0-15              10.0  12.3 

5-15               25.7  17.2 

15                7.0  4.9 

Sum                42.7  34.4 

Percentage  of  alongshore  current  episodes  *  77.1 

Percentage  of  miscellaneous  current  events  =  22.9 
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VARIATION  OF  ALONGSHORE  MEAN  CURRENT   WITH  DISTANCE  FROM  SHORE 
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Figure  1.      (After  Bull   &  Murthy  1980) 
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figure  2.      MONTHLY  VECTOR  TIME-SERES  PIOT  PICKERNG,  (After  Bull  &  Murthv.1980) 
LAKE  ONTARIO,  JANUARY  1980       -52-  V  «iuiwV,ijw< 


STATION  IE       CURRENT  VECTORS 


lw  rm  riLTR 
1»  TO  14 


I1M*) 


STATION  IE 


CURRENT   VECTORS 


I     id/*     I     ii/i      I     u/»      I     ii/f      I     u/»      I     it/»      I     !•/»     t     Wi      I     u/t      I     i»/»      I     ta/i      I     tin      I     ti/t 


lw  pmi  rn.hu 

1*  T(  14  M   CUTWT 


VECTOR  TIME-SERIES  PLOT,  DOUGLAS  POINT,  LAKE 
HURON,    JULY  28 -AUG  22,  1974 

Figure  3, 
-53- 


STATION  043 
1.3  km  FROM  SHORE 
DEPTH  6  m 


CURRENT  VECTORS 

0  10  cm/sec. 

I 1 


OFFSHORE  (135) 


W 


1/7  2/7  3/7  4/7  5/7  6/7  7/7  8/7  9/7  10/7  11/7         12/7  13/7         14/7 


STAGNATION 


SHORT  DURATION 


L 


15/7         16/7          17/7          18/7         19/7 
EPISODES 


20/7        21/7        22/7         23/7        24/7        25/7  ;*26/7        27/7         28/7 


SHORE   PARALLEL 


VECTOR  TIME-SERIES  PLOT,  LAKEVIEW,  LAKE  ONTARIO,  JULY  1  -JULY  28,1972 

Figure  4. 


MAP   OF  LAKE    ONTARIO,  PORT   CREDIT- LAKEVIEW 


metres 


0  500        WOO 


N 


PORT 
CREDIT 


LAKEVIEW 
GENERATING 
PLANT  Si: 


MARKER    FLAG 


DYE    RAFT 


CROSS    SECTION  1 


Figure  5.     Schematic  of  a  Continuous  Plume  Experiment 
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Figure  6.      crossplume  concentration  distribution -relative 
(After  Murthy  &  Miners,   1978) 
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Figure    7.        cfiossplume  concentration  distribution  -  absolute 

(After  Murthy  &  Miners,   1978) 
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Figure  9.     Schematic  of  a  laterally  diffusing  effluent  field  in 
a  steady  current. 
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Figure  9A.     Schematic  of  Plume  Dispersal  with  Different 
Diffusion  Mechanisms.     (After  Brooks,  1960) 
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Figure  10.     Coordinate  System  for  Gaussian 
Plume  Calculations. 
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Figure  11.  Schematic  of  the  Farfield  Dilution 
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FIGURE  14  DILUTION  NO  INTERFERENCE  9„  =  (7 
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FIGURE  1  5  JET  WIDTH  NO  INTERFERENCE  0.  =  0° 
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FIGURE  16  DILUTION  NO  INTERFERENCE   90  =  2CT 
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FIGURE  17  JET  WIDTH  NO  INTERFERENCE  0D  =  20T 
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FIGURE  18  DILUTION  FACTORS  FOR  AN  ALTERNATING  DIFFUSER  (STRONG  FAR-FIELD 
FRICTION  EFFECTS) 
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FIGURE  1 9  DILUTION  FACTORS  FOR  A  STAGED  DIFUSER 
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CURRENT   DILUTION 

Experimental  measurements  of  minimum  surface  dilution  of 
A  finite  line  source  of  buoyancy  flux  in  a  current. 
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Figure   20.    Current  Dilution   .(After  Roberts   1977) 
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Figure  21.  Schematic  of  the  Average  Concentration  Field 
obtained  when  Combining  Plumes  with  Different 
Orientations . 
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Figure  22.  The  average  concentration  field  of  a  single 
alongshore  current  episode. 
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Figure  23.  The  average  concentration  field  of  shore  parallel 
and  anti parallel  current  episodes. 
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Figure   23 (cont'd.) 
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Figure  24  :  Lakeview  WPCP  Outfall  Diffuser 
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A  COASTAL   DISPERSION  MODEL 
FOR  EFFLUENT  PLUMES 


APPENDIX 

Contours  of  the  mean  concentration  field  and  the  extent  (minimum, 
mean,  and  maximum)  of  the  100  concentration  contour  (when  the 
effluent  concentration  of  a  conservative  substance  was  10000)  for 
different  current  flow  regimes  are  presented  in  Figures  1.01  to  1.20. 
Monthly  summaries  and  stick  vector  plots  of  the  current  meter  data 
(including  water  temperature)   are  presented  in  Figures  1.21  to  1.31. 
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MOORING! 
START  DATE  I 
STOP   DATE i 


70-0OC-0O6 
70/  6/  8 
70/10/13 


CURRENT 

—  0-   B   CM/SEC 

—  E-1G  CH/SEC 
=>       *1S  CM/SEC 


FLOW  TOWARD 


DEPTHi  6M 


AVERAGE  CONCENTRATION  FIELD 

(BASED  ON  CURRENT  SPEED-DIRECTION  HISTOGRAM] 

CONSERVED  (SUMNER).  MEAN  LRKEVIEU  FIELD- 
SPEED  0-5  CM/SEC.   SHORE  PARALLEL . 


CM  MOORING: 
CM  OEPTHl 


70-00C-006 
6.0  H. 


START  ORTEi  70/  6/  B 
STOP  ORTEi  70/10/13 


DIFFUSIVITYi  IZOO.O  CRnZ/SEC 

DISCHARGE  RATE  i  Z100.0  L/SEC. 

DISCHAROC  SPEEDi  48-0  CM/SEC 

DENSITY  ANOMALY i  1 .5  PPT. 


SOURCE  DILUTION! 
CURRENT  SPEED ■ 
DIFFUSION  LAYERl 

OIFFUSER  DEPTHi 
OIFFUSER  WlOTHi 
OIFFUSER  ANOLEi 
SHORELINE  ANGLE i 
DIFFUSER-SHOREt 
DISTANCE-SHORE  I 


THE  CONCENTRATIONS  ARE  SCALED  BY  10--4 

THE  DIFFUSIVITY  INCREASES  A8  THE  4/3  POWER  OF 

THE  SCALE  OF  THE  DIFFUSION  FIELD. 


13.0 
4.2 
3.3 

ITD  1.0) 
CH/SEC. 
M  THICK. 

9.S 

200.0 

135-0 

225.0 

90.0 

1.3 

M. 

M. 

DEGREES. 

DEGREES. 

DEGREES. 
KM. 

OUTER  RADIUSi 
ANGULAR  INCREMENT! 
RAOIAL   INCREMENT l 


S.O   KM. 
6.0  DEGREES. 
SDO.O   METRES. 


CURRENT  HISTOGRAM l   3Z4.0   EVENTS 
SECTOR   INCREMENT!    S.O  DEGREES. 
SPEED   INCREMENT!     1.0  CH/SEC. 


Figure  1.01       Contours  of  the  Mean  Concentration  Field 
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Figure  1.01      Contours  of  the  Mean  Concentration  Field 
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CURRENT  RDVECTION  DIAGRAM   * 

OFF      1 1 36 1 


II 


ANT  I    I    «S) 


(3^ 


PflRfi    .225) 


ONTO   (316) 


affMSTMsraussr- 

STUB  WRE  —  "  ■ 


III 


CURRENT  HISTOGRAM 

(32*  EVENTS.   6.0'  SECT0R6) 


*  *  » 


HOORINO* 
START  DATE. 
STOP  DATE • 


70-00C-006 

70/  •/  • 
70/10/13 


CURRENT 

0-6  CM/SEC 

m—    6- IE   Ctl/SEC 

c=      215  Cn/SEC 


FLOW   TOMBRO 


„»      1  ONSHORE 
K  11361 


PARALLEL 
(226) 


OEPTH.    6« 


III 


ouFRflr,F  CONCENTRATION  FIELD 

flyB«EOH0;CUSsPEEoioiRECT.(mHI8TO«.RH. 

CONSERVED   (SUH-IER..   0.01   CWJOUO  «««. 
LBKEVIEN    H1N.    HERN.    MAX .    10-5    CM/S.      I 


III 


CM  MDORINO.  7D-0QC-0O6 
en  depth:    8-o  n. 


START  OATEi  70/  6/  8 
STOP  OATEi  70/10/13 


OIFFUSIVITY.   o  1800.0  CM..2/SEC . 

DISCHARGE  RATE.  2100.0  ^/SEC. 

D15CMRR0E  SPEEDi  «-°  "t 

DENSITY  ANOMALY.  I •«  rr\ ■ 

SOURCE  DILUTION.     tl-O   (TOl-Ol 
CURRENT  SPEED.        «•*  "/BEC. 
DIFFUSION  LATER.      3.9  "  ™1C*' 

DIFFUSER  DEPTH. 
01FFUSER  WIDTH. 
DIFFUSER  BNOtE: 
SHORELINE  BNOLE. 
D1FFUSER-SH0RE. 
DISTANCE-SHORE ' 


9.6 

200.0 

135.0 

226.0 

90.0 

1.3 


ri- 
ll. 

OEOREES. 
DEGREES. 

OEOREES- 
KM. 


T«l  SCALE  OF  THE  DIFFUSION  F!ELD. 

iNCRtneNTi  soo.o  metres- 


RADIAL 

CURRENT  H1STQ0RBH. 
SECTOR  INCREMENT. 
SPIED   INCREMENT. 


324.0  EVENTS 
5.0  DECREES. 
1.0  CM/SEC 


Figure  1.02        Extent  of  the  100  Concentration  Contour 
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Figure  1.02        Extent  of  the  100  Concentration  Contour 
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»»»   CURRENT  flDVECTION  DIAGRAM   •*» 

OFF   (136 1 


in   CURRENT  HISTOGRAM   «■■ 

(420  EVENTS.   6.B*  SECTORS) 


BUT  I  (  46! 


MOOR  I NO ■     70-00C-0O6 
START  DATEl   70/  7/22 
STOP  DBTEt  70/10/24 


PARA  (226) 


ONTO  1316) 


THE  IRAEQULAA  CONTOUR  CIRCUMSCRIBES  THE  ""HATED 
HAXinun  ADVECTION  OF  THE  SECTOR  CURRENT  SUMS. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  10.0  KM 
RADIUS  CIRCLE  OF  THE  MAIN  PLOT. 


CURRENT 

0-6  CH/6EC 

—  6-16  CH/SEC 

>==  us  cn/sec 


FLOH  TOWARD 


PARALLEL 
225) 

*X0KSM0RE 
( 1351 


I6X 


DEPTH  ■  OH 


,„   AVERAGE  CONCENTRATION  FIELD 

IBRSED  ON  CURRENT  SPEED-DIRECTION  HISTOORRHI 


■  KM 


CONSERVED  (SUNHER) 
SPEED  6-16  CM/SEC 


MEAN  LflKEVIEH  FIELD. 
SHORE  PARALLEL. 


CM  HOORINOi  70-00C-006 
CH  OEPTHi     6-0  H. 


STAltr  OATEi  70/  7/ZZ 
STOP  OBTEl  70/10/24 


OJFFUSIVlTYi  1200.0  CN«2/SEC 

01SCHBROE  RBTEi  2100.0  L/SEC . 

01SCHRROC  SPEEOi  48.0  CH/SEC. 

DENSITY  ANOHRLTi  1-6  PPT. 


SOURCE  DILUTION! 
CURRENT  SPEEOi 
DIFFUSION  LBYERi 

OIFFUSER  OEPTHi 
DIFFUSER  MIDTMi 
OIFFUSER  ANOLEi 
SHORELINE  BNOLEi 
OlFFUSER-SHOREi 
DISTANCE-SHORE  I 


13.0   (TO  1.0) 
6. 9  CH/SEC. 
2.1  M  THICK. 


9.6 

200.0 

136.0 

226.0 

90.0 

1.3 


H. 

n. 

OEOREES. 
DECREES. 
OEOREES. 

KM. 


.  THE  CONCENTRATIONS  RRE  SCALED  BT  10--4 
.  THE  01FFUS1VITY  INCREASES   AS  THE  4/3  POWER  OF 
THE  SCALE  OF  THE  DIFFUSION  FIELD. 

OUTER  RADIUS i        10.0  ««■ 
BNOULAR  INCREMENT i    6.0  OEOREES. 
RROIBL   INCREHENTl  1000.0   METRES. 


CURRENT  HISTOGRAM: 
SECTOR  INCREMENT i 
SPEED    INCREHENTl 


420.0  EVENTS 
S.O  DECREES- 

1.0     CH/SEC. 


Figure  1.03 


Contours  of  the  Mean  Concentration  Field 
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Figure  1.03  Contours  of  the  Mean  Concentration  Field 


-83- 


...   CURRENT  ADVECTION  DIAGRAM   km 

OFF   1136) 


ANT  I  I  45) 


*kk   CURRENT  HISTOGRAM   ««* 

[4Z0  EVENTS.  S.O'  SECTORS! 


STRKT  DRTEl 
STOP  OATEt 


PftRfl  (225) 


ONTO  13)6! 


THE  IRREOULAH  CONTOUR  CIRCUMSCRIBES  THE  ESTIMATED 
MAXIMUM  ADVECTION  OF  THE  SECTOR  CURRENT  SUMS. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  10.00  KM 
RAOIUS  CIRCLE  OF  THE  MAIN  PLOT. 


CURRENT 

—  0-  B  CM/8EC 

—  6-16  CM/SEC 
=a   *16  CH/SEC 


FLOM  TOHARO 


DEPTHi    6H 


III 


AVERAGE  CONCENTRATION  FIELD 

IBASED  ON   CURRENT   SPEEO-DIRECT ION  HISTOORAH) 


■  KM 


CONSERVED  (SUMMER) 
LAKEV1EM  HIN.  MEAN 


0-01  CONTOUR  RANBE. 

nAx.  (6- is  cn/s.  P) 


CH  MOORING-)  70-OOC-006 
CH  DEPTHi      6.0   M. 


START  DATE i  70/  7/22 
STOP  DATE  I  70/10/24 


DIFFUSIVITTi  12D0.0  CM.-2/SEC 

DISCHARGE  RATE i  2100.0  L/SEC . 

DISCHARGE  SPEEDi  4B.0  CH/SEC. 

DENSITr  ANOMALY)  1.6  PPT- 


60URCE  DILUTION. 
CURRENT  SPEEDi 
DIFFUSION  LAYER) 

DIFFUSER  DEPTHi 
DIFFUSE*  WIDTH) 
DIFFUSER  ANOLE< 
SHORELINE  ANOLE) 
DIFFUSER -SHORE i 
DISTANCE -SHORE' 


13-0   (TO  1.0) 
6.9  CH/SEC 
2.7  H  THICK. 


9.6 

200-0 

136. 0 

226. 0 

90.0 

1.3 


H. 

H. 

DEOREES- 

DEOREES. 

DE0REE6. 

KM. 


THE  CONCENTRATIONS  ARE  SCALED  BY  10"4 

THE  01FFUS1VITY  INCREASES  AS  THE  4/3  POMER  OF 

THE  SCALE  OF  THE  DIFFUSION  FIELD. 


OUTER  RADIUS) 
ANOULAR  INCREMENT  i 
RADIAL   INCREMENT) 

CURRENT  HlSTOORAni 
SECTOR  INCREMENT) 
SPEED    INCREMENT) 


10.00 
6-0 

1000. 0 


KH. 

DEGREES. 

METRES. 


420.0   EVENTS 
S.O  OEDREES. 
1.0  CH/SEC. 


Figure  1.04 


Extent  of  the  100  Concentration  Contour 
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Figure  1.04  Extent  of  the  100  Concentration  Contour 
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in   CURRENT  RDVECTION  DIAGRAM   ■■» 

OFF      (13St 


RNTI    (    45) 


O 


im      CURRENT    HISTOGRAM      ■  ■■ 

(342  EVENTS.   E.Oa  SECTORS) 


MOOR  I  MOt     70-000006 
START  DRTCt  70/  6/26 
•  TOP  DRTEr   70/10/31 


PARR  IZZS I 


ONTO  (316) 

THE  IRREOULAR  CONTOUR  CIRCUMSCRIBES  THE  ESTIMATED 
MAXIMUM  RDVECTION  OF  THE  SECTOR  CURRENT  SUNS. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  5 .0  KM 
RADIUS  CIRCLE  OF  THE  HAIN  PLOT. 


CURRENT 

0-5  CH/SEC 

6-16  CH/SEC 

<=   *16  CH/SEC 


FLOW  TOHRRO 


OEPTHi  6H 


III 


AVERAGE  CONCENTRATION  FIELD 

(BASED  ON  CURRENT  SPEED-DIRECTION  HISTOORAH) 


MX* 


CONSERVED  (SUtlHER) 
SPEEO  0-5  CH/SEC. 


HERN  LAKEVIEM  FIELD. 
SHORE  ANTIPRRRLLEL. 


Cfl  HOORINOt  70-00C-D06 
CM  OEPTHi     6.0  n. 


START  DATE  I  70/  5/26 
STOP  DATEi  70/10/31 


OlFFUSIVITTi  1200.0  CH«»2/SEC . 

01SCHAROE  RATEi  2100.0  L/SEC. 

DISCHRROE  SPEED  i  46.0  CH/SEC 

DENB1TV  ANOMALTi  1.6  PPT. 


SOURCE  DILUTION! 
CURRENT  SPEEO  ■ 
DIFFUSION  LAYER  I 

DIFFU6ER  OEPTHi 
DIFFUSER  NIOTHi 
D1FFUSER  ANGLE  I 
SHORELINE  ANOLEi 
01 FF USER-SHORE i 
DISTANCE-SHORE  i 


13.0   (TO  1.0) 

4.1   CH/SEC. 
4-1   H  THICK- 


9-6 

200.0 

136-0 

226. 0 

SO.O 

1.3 


H. 

n. 

DEOREES. 
DEOREES. 
DEOREES. 
KM. 


THE  CONCENTRATIONS  ARE  SCALEO  BT  I0»4 

THE  DIFFU51V1TT  INCREASES   AS  THE  4/3  POWER  OF 

THE  SCRLE  OF  THE  01FFUSIDN  FIELD- 


OUTER  RAOIUSi 
AMHJLAR  INCREHENTl 
RADIAL   INCREHENTl 

CURRENT  HISTOORAH i 
SECTOR  INCREHENTl 
SPEEO    INCREHENTl 


6-0 

s.o 

600.0 


KH. 

OEOREES. 

nETRES. 


342.0  EVENTS 
S-0  DEOREES. 
1.0  CH/SEC. 


Figure  1.05 


Contours  of  the  Mean  Concentration  Field 
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Figure  1.05  Contours  of  the  Mean  Concentration  Field 
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.»«   CURRENT  RDVECTION  DIRGRfiM   «*« 

OFF   1138) 


PARA  I  226) 


ANT.  ,  46.  <^T^?  1 

ONTO  13161 

THE  IAREOULAR  CONTOUR  CIRCUMSCRIBES  THE  ESTIMATED 
MAXIMUM  ADVECTtON  OF  THE  SECTOR  CURRENT  SUMS. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  6.0  KM 
RADIUS  CIRCLE  OF  THE  MAIN  PLOT. 


*■«   CURRENT  HISTOGRRM   ■  ■* 

(S42   EVENTS.      S.0«    SECTORS) 


MOORINQt  7Q-00C-0OB 

START   DATE  I      70/   6/28 
STOP      DATE i     70/10/31 


CURRENT 

0-6  CM/SEC 

6-16  CM/SEC 

■a      i!6  CH/SEC 


FLOW   TOMARO 


DEPTH  I    6M 


*** 


AVERAGE  CONCENTRATION  FIELD 

(BASED  ON  CURRENT  SPEED-DIRECTION  HISTOGRAM) 

CONSERVED  I  SUMMER ) 
LAKEVIEW  HIN.  MEAN 


«»■ 


0.0!  CONTOUR  RANOE. 
MAX.  10-5  CM/S.   AP) 


CM  MOOR  I  NO 

CM  DEPTH i 


7Q-00C-006 

e.o  n. 


START  DATE i  70/  6/2B 
STOP  DATE  i  70/10/31 


OIFFUSIVITTi  1200.0  CMm2/SEC. 

OISCHAROC  RATE  ■  2100.0  L/SEC 

DISCMAROE  SPEED'  40  .0  CM/SEC 

DENSITY  ANOMALTi  1.6  PPT. 


SOURCE  DILUTION' 
CURRENT  SPEED i 
DIFFU8I0N  LATERi 

DIFFUSER  DEPTH i 
OIFFUSER  HIDTHi 
DIFFUSER  ANSLEi 
SHORELINE  ANGLE i 
OlFFUSER-SHOREl 
DISTANCE-SHORE  I 


13.0   (TO  1.0) 
4.1   CM/SEC 

4.1   M  THICK. 


8.6 

200.  D 

136.0 

226.0 

90.0 

1.3 


n. 

M. 

DEOREES. 

OEOREES. 

DEOREES. 
KM. 


THE  CONCENTRATIONS  ARE  SCALEO  BT  10.-4 

THE  OIFFUSIVITT  INCREASES  AS  THE  4/3  POWER  OF 

THE  SCALE  OF  THE  DIFFUSION  FIELD. 

OUTER  RAOIUBi         6.0  KM. 
ANDULAR  INCREMENT t    6.0  OEOREES. 
RADIAL   INCREMENT!   600.0  METRES. 

CURRENT  HISTOORAMi   342. 0  EVENTS 
SECTOR   INCREMENTi    8.0  OEOREES. 
SPEED    INCREMENTi     1.0   CM/SEC. 


Figure  1,06    Extent  of  the  100  Concentration  Contour 
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Figure  1.06    Extent  of  the  100  Concentration  Contour 
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«««   CURRENT  HDVECTION  DIAGRAM   *** 

OFF  U35) 


ANT)  I  4S) 


«**   CURRENT  HISTOGRAM   *«* 

(528  EVENTS.   5.0'  SECTORS) 


MOOR  I  NO ■ 
START  DATE! 
STOP  OBTEi 


PARA  12261 


ONTO  I31S) 


THE  IRREGULAR  CONTOUR  CIRCUMSCRIBES  THE  ESTIMATED 
HRXH1UK  ADVECTION  OF  THE  SECTOR  CURRENT  SUMS. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  10-0  KM 
RADIUS  CIRCLE  DP  THE  MAIN  PLOT. 


CURRENT 

—0-5  CM/SEC 
—  6-16  CM/SEC 
=   >1B  CM/SEC 


FLOW  TONARO 


DEPTH!  6n 


OVERAGE  CONCENTRATION  FIELD   ■■■ 

(BASED  ON  CURRENT  SPEEO-OIRECTION  HISTOGRAM) 


CONSERVED  (SUMMER) 
SPEED  6- IS  CM/SEC. 


MEAN  LAKEVIEM  FIELD. 
SHORE  ANTIPARALLEL. 


CH  MOORlNOi  70-0QC-0D8 
CM  DEPTHi     60  H. 


START  DATE i  70/  5/2* 
STOP  OATEf  70/10/29 


DIFPUSIVITVi  1200.0  CH«»2/SEC. 

OISCHAROE  RATE:  2100.0  L/SEC . 

DISCHARGE  SPEEDi  48.0  CM/SEC. 

0ENS1TT  ANOHALVi  1.6  PPT. 


SOURCE  DILUTION! 
CURRENT  SPEEDi 
DIFFUSION  LAYER! 

DIFFUSER  DEPTH  ■ 
DIFFUSER  WlOTMi 
DIFFUSER  ANDLEi 
SHORELINE  ANGLE i 
DIFFUSER-SHOREi 
DISTANCE-SHORE l 


THE  CONCENTRATIONS  ARE  SCALED  BY  10a*4 

THE  QIFFUSIVITY  INCREASES   AS  THE  4/3  PONER  OF 

TME  SCALE  OF  THE  DIFFUSION  FIELD. 

OUTER  RADIUS!        10.0  Kn. 
AN'OULAR  INCREMENT  I    6.0  OEOREES . 

RADIAL   INCREMENT!  1000.0  METRES. 

CURRENT  HISTODRAMl   5ZB.D   EVENTS 
SECTOR   INCREMENT!    6.0  OEOREES. 
SPEED    INCREMENT!    1.0  CM/SEC 


13.0 

(TO   1.0) 

6-9 

cm/sec 

2.6 

M   THICK. 

9.6 

M. 

200.0 

H. 

136.0 

OEOREES. 

226.0 

OEOREES. 

90.  C 

OEOREES. 

1.3 

KM. 

Figure  1.07 


Contours  of  the  Mean  Concentration  Field 
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Figure  1.07  Contours  of  the  Mean  Concentration  Field 
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***   CURRENT  ADVECTION  DIAGRAM  *** 

OFF      1 1  35  1 


ANTI    (    46) 


*€> 


**«      CURRENT    HISTOGRAM      *** 

I52B  EVENTS.   5.Q*  SECTORS  I 


HOOD I NO ■ 
START  DATE! 
STOP  DATEi 


PARR  (226) 


ONTO  1 31E ) 

THE  IRREOULAR  CONTOUR  CIRCUMSCRIBES  THE  ESTIMATED 
MAXIMUM  ADVECTION  OF  THE  SECTOR  CURRENT  SDKS. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  10.00  KM 
RROIUS  CIRCLE  Of  THE  BAIN  PLOT. 


CURRENT 

—  0-  6  Ctl/SEC 

—  5-16  CH/SEC 

=   HE  CM/SEC 


FLOW  TOWARD 


DEPTH  i  6H 


MM  * 


AVERAGE  CONCENTRATION  FIELD 

I  BASED  ON  CURRENT  SPEEO-DIRECTION  HISTOORAH) 


CONSERVED  (SUMMER)  . 
LAKE VI EH  MIN.  MEAN, 


0-01  CONTOUR  RANOE. 
MAX.  15-15  CM/S.  API 


CM  MOORINOi  70-00C-OO6 
CM  DEPTHi     6.0  A. 


START  DATE>  70/  5/24 
STOP   DATEt  70/10/28 


DIFFUSIVITr I 

1200-0 

CM.<2/SEC 

OISCHAROE  RATEi 

210D.0 

L/SEC. 

OISCHAROE  SPEED i 

4B.D 

CM/SEC. 

OENSITT  ANOMALY! 

1.5 

PPT. 

SOURCE  DILUTION! 

13.0 

(TO  1.0) 

CURRENT  SPEED) 

6.9 

CM/6EC. 

DIFFUSION  LATER! 

2.6 

N  THICK. 

0IFFU8ER  DEPTHi 

B.5 

M. 

DIFFUSER  HIDTHt 

200.0 

n. 

DIFFUSER  AKJLEi 

136. 0 

OEOREES. 

SHORELINE  ANOLE ■ 

226-0 

OEOREES. 

DIFFUSER-SHOREl 

90.0 

OEOREES. 

OlSTANCE-SHOREi 

1.3 

MM. 

THE  CONCENTRATIONS  ARE  SCALED  Br  10«4 

THE  0IFFUSIV1TT  INCREASES   AS  THE  4/3  POWER  OF 

THE  SCALE  OF  THE  DIFFU6I0N  FIELD. 


OUTER  RADIUS! 
ANOULAR  INCREMENT! 
RAOIAL   INCREMENT i 


10.00  KM. 

6-0   OEOREES. 
1000.0   METRES. 


CURRENT  HISTOGRAM i   628.0  EVENTS 
SECTOR   INCREMENT!     S.O   OEOAEES. 
8PEE0    INCREMENT!     1.0  CM/SEC. 


Figure  1.08 


Extent  of  the  100  Concentration  Contour 
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Figure  1.08    Extent  of  the  100  Concentration  Contour 


-93- 


kxm   CURRENT  ADVECT10N  DIAGRAM   *** 

OFF   (13SI 


ANT  I  (  45  J 


ONTO  1315) 


***   CURRENT  HISTOGRAM   *** 
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Contours  of  the  Mean  Concentration  Field 
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Figure  1.09  Contours  of  the  Mean  Concentration  Field 
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Figure  1.10  Extent  of  the  100  Concentration  Contour 


■97- 


...   CURRENT  ADVECTION  DIAGRAM   «»« 

OFF   1136) 


ANT  I  I  46) 


»*«   CURRENT  HISTOGRAM   «■■ 

(604  EVENTS.   6-0'  SECTORS) 


tlOORINDi     SEVERAL 
START  DATE i   70/11/  6 
STOP  DATE  I   71/4/1 


PARA  1226) 


ONTO  (315) 
THE  IRREOULAR  CONTOUR  CIRCUMSCRIBES  THE  ESTIMATED 

maximum  aoveciion  of  the  sector  current  suns. 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  10-0  M 
RADIUS  CIRCLE  OF  THE  MAIN  PLOT. 


CURRENT 

0-6  CM/SEC 

—  E-IS  CH/SEC 
<=»   116  CM/SEC 


FLOW   TOWARD 


OEPTMl    en 


111 


AVERAGE  CONCENTRATION  FIELD 

(BASED  ON  CURRENT  SPEED-DIRECTION  HISTOGRAM  I 


CONSERVED  (K INTER  I 
SPEED  6-16  CM/SEC. 


MEAN  LAKEVIEN  FIELD. 
SHORE  PARALLEL. 


CM  MOORINDi  SEVERAL 
CM  DEPTMi     6.0  H. 


START  DATEi  70/11/  6 
STOP  OATEi  71/  4/  S 


DIFFUSIVITTl 

1200. 0 

CH»»2/SEC 

DISCHARGE  RATE: 

2100.0 

L/SEC 

OISCHARDE  SPEED i 

48.0 

CH/SEC. 

OEMS  ITT  RN0I1AIV< 

1.6 

PPT. 

SOURCE  DILUTION! 

13.0 

(TO  1-0) 

CURRENT  SPEEOi 

8.3 

CM/SEC. 

01FFUSION  LATER! 

2.4 

H  THICK. 

DIFFUSFR  OEPTHi 

9.6 

M. 

DIFFUSE**  MlOTHl 

ZOO.  0 

H- 

DIFFUSER  ANOLE 1 

135.0 

DECREES. 

SHORELINE  ANOLE i 

225.0 

GEOREES. 

DIFFUSER-SMOREl 

90-0 

OEOREES. 

DISTANCE-SHORE  I 

1.3 

KM. 

THE  CONCENTRATIONS  ARE  SCALED  BT  10»»4 

THE  OIFFUSIVITT  INCREASES  AS  THE  4/3  POWER  OF 

THE  SCALE  OF  THE  DIFFUSION  FIELD. 

OUTER  RADIUS!        10.0  KM. 

HNOULAR  INCREMENT i    5.0  OEOREES. 

RADIAL   INCREMENT!  1000.0  HETRES. 

CURRENT  HISTOORAHt   604.0  EVENTS 

SECTOR   INCREMENT!    5.0  DECREES. 

SPEED   INCREMENT!     1.0  CM/SEC. 


Figure  1.11 


Contours  of  the  Mean  Concentration  Field 


-98- 


Figure  1.11  Contours  of  the  Mean  Concentration  Field 
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Figure  1.12         Extent  of  the  100  Concentration  Contour 
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Figure  1.13  Contours  of  the  Mean  Concentration  Field 
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Figure  1.13  Contours  of  the  Mean  Concentration  Field 
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Figure  1.14  Extent  of  the  100  Concentration  Contour 
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Contours  of  the  Mean  Concentration  Field 
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Figure  1.15  Contours  of  the  Mean  Concentration  Field 
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Extent  of  the  100  Concentration  Contour 
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Figure  1.17  Contours  of  the  Mean  Concentration  Field 
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Figure  1.18    Extent  of  the  100  Concentration  Contour 
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Figure  1.19  Contours  of  the  Mean  Concentration  Field 


■115- 


...   CURRENT  ADVECTION  DIAGRAM   «** 

OFF   (1361 


ANTI  I  *S  1 


0 

ONTO  (315) 


...   CURRENT  HISTOGRAM   *«* 

(95  EVEHTS.   6.0«  SECTORS] 


MOOR  I  NO ■ 
STMT  ORTEi 
STOP  OATEt 


PARA  (22St 


CURRENT 

O-  G  CN/SEC 

—  B-1S  CH/SEC 

==  »i5  cn/SEC 


THE  IRHEOULHR  CONTOUR  CIRCUMSCRIBE6  THE  ESTIMATED 
MAXIMUM  ADVECTION  OF  THE  SECTOR  CURRENT  SUNS- 

THE  ABOVE  CIRCLE  REPRESENTS  THE  OUTER  20.0  KM 
RADIUS  CIRCLE  OF  THE  MAIN  fLOT. 


FLOW  TOWARD 


OErTHi  6M 


...   AVERAGE  CONCENTRATION  FIELD 

(BASED  ON  CURRENT  SPEED-OIRECTION  MI8T00AAMI 

CONSERVED  (HINTER).  0.01  CONTOUR  RANOE . 
LAKEVIEN  NIN.  MEAN.  HAX .  I>  IB  CH/S,  API 

CM  MOORINO.  71-00C-008     START  DATE.  71/  2/21 
CM  DEPTH.     6.0  N.      STOP  OATEi  71/  3/11 

DIFFUSIVIRi  1200.0  CM«2/SEC. 

OISCHAROE  RATE  I  2100.0  L/SEC 

DISCMAROE  SPEED!  48-0  CM/SEC . 

DENSITY  ANOMALY «  1.6  PPT. 

SOURCE  DILUTION! 
CURRENT  SPEEOi 
DIFFUSION  LATER! 

DIFFUSEH  DEPTH! 
OIFfUSER  HIDTMt 
DIFFUSER  ANOLEi 
SHORELINE  ANOLEi 
OIFFUSER-SMOREl 
DISTANCE-SHORE  I 

.  THE  CONCENTRATIONS  ARE  SCALED  8T  »0"« 
■  THE  DIFFUSIVITY  INCREASES  AS  THE  4/3  POMER  OF 
THE  SCALE  OF  THE  01FFUSI0N  FIELD. 

OUTER  RAOIUSi        20-0  "H. 
ANOULAA  INCREMENT!    B.D  DEOREES. 
RACIAL   INCREMENT!  2000.0   METRES. 

CURRENT  MISTOORAHI    9S.0   EVENTS 
SECTOR   INCREMENT!     S.O   DECREES. 
SPEEO    INCREMENT!     1.0   CM/SEC. 


... 


13.0 

1*.* 

1.4 

(TO    1-0! 
CM/SEC 
H   THICK. 

9.5 

n. 

200.0 
136-0 

226.0 
90-0 

Hi 

DEOREES. 
DEOREES. 
DEOREES. 

1.3 

KM. 

Figure  1.20  Extent  of  the  100  Concentration  Contour 


■116- 


10 


,,.   !*•   i6* 


iv"     »■ 


Figure  1.20    Extent  of  the  100  Concentration  Contour 
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PARALLEL 
(225) 


ONSHORE 
(13B) 


DEPTH «    6   H 
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Figure   1.31 


NOTES 
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